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Foreword 


I thank Springer-Verlag for keeping my Introduction to Diophantine 
Approximations in print. This second edition is unchanged from the first, 
except for the addition of two papers, written in collaboration with 
W. Adams and H. Trotter, giving computational information for the 
behavior of certain algebraic and classical transcendental numbers with 
respect to approximation by rational numbers and their continued frac- 
tions. I thank both of them for their agreement to let me reproduce 
these papers, which expand and illustrate the general theory in computa- 
tional directions. 

The classical numbers, as I described them in 1965, are those which 
can be obtained by starting with the rational numbers, and performing 
the following operations: 


~ Take the algebraic closure, thus obtaining a field F. 

— Take a classical, suitably normalized transcendental function (elliptic, 
hypergeometric, Bessel, exponential, logarithm, etc.), or jazzed up ver- 
sions, coming from normalized transcendental parametrizations of alge- 
braic varieties, take values of such functions with argument in F, and 
adjoint them to F. 

— Iterate these two operations inductively. 


Questions arise as to the properties of the numbers so obtained (a 
denumerable set), from the point of view of diophantine approximations. 
The present book may be viewed as providing the simplest examples at 
the most elementary level, using only the most elementary language of 
mathematics. 


New Haven, 1995 SERGE LANG 


if 


es 


Foreword to the First Edition 


The quantitative aspects of the theory of diophantine approximations are, 
at the moment, still not very far from where Euler and Lagrange left 
them. Very recent work seems to have opened some fruitful lines of 
research, and in this book we shall illustrate by significant special exam- 
ples three aspects from the theory of diophantine approximations. 

First, the formal relationships which exist between various counting 
processes and functions entering in the theory. These essentially occur in 
Chapters I, II, II. 

Second, the determination of these functions for numbers which are 
given as classical numbers, in a concrete fashion. Chapters IV and V 
give examples of this. 

Third, we have mentioned certain asymptotic estimates holding almost 
everywhere (e.g. the Khintchine theorems and the Leveque—Erdés— 
Schmidt theorems). Such results are useful since they suggest roughly 
what may be considered “pathological” numbers, and also the range of 
magnitude of similar estimates for the classical numbers. However, as 
one sees from the quadratic numbers (which are of constant type), and 
the Adams result for e, each special number may exhibit its own particu- 
lar behavior in the more subtle range of approximation. To determine 
this behavior for the classical numbers is perhaps the most fascinating 
part of the theory of diophantine approximations. 

There exist other aspects, for instance the connection with transcen- 
dental numbers, but these have been left out completely since the style of 
the results known in this direction is at present so different from the style 
of the results which we have emphasized here. 

I have avoided including partial results whose statements seemed to 
me too remote from expected best possible statements. Every chapter 


Viii FOREWORD TO THE FIRST EDITION 


should be viewed as working out a special case of a much broader 
general theory, as yet unknown. Indications for this are given through- 
out the book, together with references to current publications. 

It is unusual to find a mathematical theory which is in a state as 
primitive and naive as the present one, and there is of course some 
delight in catching it in that state. In fact, this book may be used for a 
course in number theory, addressed to undergraduates, who will thus be 
put in contact with interesting but accessible problems on the ground 
floor of mathematics. If, however, like Rip van Winkle, I should awake 
from slumber in twenty years, my greatest hope would be that the theory 
by then had acquired the broad coherence which it deserves. 


Berkeley, 1966 SERGE LANG 
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CHAPTER | 


General Formalism 


I, §1. RATIONAL CONTINUED FRACTIONS 


We are interested in the following problem. Given an irrational number 
a, determine all solutions of the inequality 


(1) We Hee 
q 


with integers q, p, or more generally, of an inequality 


(2) |qa — p| < w(q), 


where y is some positive decreasing function of a real variable. If € is a 
real number, we may write |¢|| for the distance between & and the 
nearest integer. Then |qa« — p| = ||qa|| whenever it is sufficiently small, 
and we may rewrite the fundamental inequality (1) in the form ||qa|| < 1/q. 
In investigating ||qa|], we are therefore interested only in the residue class 
of qx modulo Z, where Z is the additive group of integers. The factor 
group R/Z (R = real numbers) is sometimes called the circle group, since 
it is isomorphic to the group of complex numbers of absolute value 1 
under the mapping x +> e?”*. We may think of || || as a metric on the 
circle R/Z. 

The inequality (1) plays a fundamental role, and it turns out that one 
can describe most of its solutions by a rational process. In this chapter, 
we shall describe this process. The more general inequality (2) will be 
considered in the next chapter. 
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Unless otherwise specified, the results of this chapter are due to Euler 
and Lagrange. (For specific references, cf. Perron [21], which contains an 
extensive bibliography of the older literature, and is an excellent reference.) 

We start by considering independent variables ay, a,, a ,.... We shall 
define inductively pairs of polynomials 


Pr = Pr(@o, esx sd.) and daa An(Ao» co sta); 


starting with po = dy and qg = 1. The quotient p,/q, will be written 


Pp. 
= [a, Beas || 
n 


Suppose that we have defined p,, gq, with k <n, and n21. We shall 
use the abbreviation 


Pk = Pr(G1,-+-5,) and = gy = Gy (1, ..- Qn). 
We define inductively 
(3) Py =O6Pn 21 + a1 and Qn = Dn-1+ 


This implies that 


Pn 1 
(4) [ao,-.-,4,] =— =a FSS SS 
Qn [a,, Seed | 
which, written in full, is equal to 
Ay + : 
i @la sp : 
a, + : 
a a; + 
QAn-1 Lee 
a 


n 


The sequence of fractions {p,/q,} is called a continued fraction. When we 
substitute numbers for ay, a,,... such that q, does not vanish, we obtain 
a sequence of numbers, which is still called a continued fraction. 


Theorem 1. We have for n= 2, 


Pas GnPn-1 ae Pn-2> 
an == an@n-1 SP Qn-2- 
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Proof. For n= 2, the assertion is verified directly. Assume n > 2, and 
assume inductively that 
Pr-1 = GnPn-2 + Pn—s> 
Qn—1 = AnGn—2 + In—3: 


Using (3), we find 


Pn = 49(AnPn—2 + Pn—3) + GnGn—-2 + In-3 
= An(4oPp—2 + In—-2) + AoPr—3 + Gp-3 


= GyPn-1 at Pn-2 
and 


Qn = An,Pn—2 ap Pn-3 = AnAn-1 ak Qn-2> 
thereby proving our theorem. 


For convenience, we define p_,=1 and g_,=0. Then Theorem 1 
remains valid when n = 1. 

In applications, we shall be interested in the nature of the values of p, 
and gq, when we substitute real numbers, for ay, a,,.... We shall always 
take a,,a,... to be >0. In that case, we see inductively that g, > 0 for 
n= 1, and hence that the fraction p,/q, has meaning as a real number. 
In particular, we obtain a corollary for Theorem 1. 


Corollary 1. Let a; be real numbers >0 for 1 Sign. For1Skn, 


let 
pe Weel 
Then 
la er NiGe...-52p=15 he) 
— Pu-itk + Pr-2 
Q-1% + W-2 
Corollary 2. Let do,...,a, and bo,...,b, be real numbers such that 


a;=1 for i21 and also b,21 for i21. Assume that a;,b; are inte- 
gers forO<sjsn-1. If 


[ao; Wee 4 = [bo ee 
then a; = b; for i 2 0. 
Proof. Jet r, = [4;,-.-,4, |, so that 


1 


r,; =a, + —_—_.. 
; ere aa 
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Then r, = 1, and similarly, s, = [b,,...,b,] 2 1. By hypothesis, 


lf 7, = 1, then 


is an integer, and hence 1/s, is also an integer. Hence s; = 1 and ay = Do. 
If r; > 1, then 


ag + — 
0 ry 


is not an integer, and hence s, > 1 also. But then a) = by because both 
Qo, bo are the greatest integers <[ap,...,a,]. Thus in all cases, ag = bo 
and r, = s,. We can now conclude the proof by induction. 


Corollary 2 gives us a uniqueness for the continued fraction formed 
with real numbers under the hypotheses of this corollary. It will be 
applicable in the next section. 

For the next theorem, we return to indeterminates. 


Theorem 2. For n = 0 we have 


QnPn-i — PnAn-1 = = We 


Proof. We have 
qoP-1 — Pog-1 = 1. 


In general, multiply the first expression in Theorem 1 by q,_,, multiply 
the second by p,_,, and subtract the first from the second. We obtain 


QnPn—-1 — PnQn-1 = sai? Pes eee — Pn-1 Qn-2)- 
This proves the theorem, because it shows that when n changes by one 


unit, the expression on the left of the inequality in the theorem changes 
by a minus sign. 


Corollary 1. For n= 1 we have 
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We shall be interested in the values taken by p, and q, when dy, ay,... 
are integers. We shall assume throughout that when we substitute such 
integers, then a,,a,... are always > 0. 


Corollary 2. If a,,a,... are positive integers, then p, and q, are 
relatively prime, and 


is ae se 
forms a strictly increasing sequence of integers. 


Corollary 3. Let « denote the rational function 


“= [ay, Teen eee b 
Then 


(— 1yrs 


Gre Cian Dr a - 
a i. Qn+29nt1 + Qn 


Proof. Replace n by n+ 2 in the theorem, and divide by q,,,. Note 
that by definition, p,.>/¢,+. = % and 9,42 = G)429n+1 + 4,- The relation 
of our corollary then drops out. 

Theorem 3. For n= 1 we have 

QnPn—2 — PnQn-2 = (= 1 ta 


Proof. We multiply the first expression in Theorem 1 by q,_,, multi- 
ply the second by p,_, and subtract the first from the second. We 
obtain, using Theorem 2, 


QnPn—2 — PaQn-2 = An(Qn—1 Pn-2 ba jet: Be, = (—1}'a., 
as was to be shown. 


Corollary 1. For n 2 2 we have 


Pn-2 ars Pn ns = 1 ter 
Qn-2 Qn GnIn-2 


Corollary 2. If a,,a,... are positive numbers, then the sequence p,/q, 
for even n is strictly increasing, and for odd n, it is strictly decreasing. 


Corollary 3. Let « denote the rational function 
% = [do,.--,4y+2]. 


(—1)"a,42 


Gn4+29n+1 aa In 


Then 


Qn® — Pa = 
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Proof. Replace n by n+2 in Theorem 3, and divide by q,+2. Note 
that by definition, Pr2/dn+2 = % ANd Ga2 = 4n+29n+1 + 4n- The relation 
of our corollary then drops out. 


Theorem 4. For n= 1 we have 


Qn 
Qn-1 


= gna cy |. 


Proof. For n=1, the assertion is clear. Assume n> 1. Suppose 
inductively that we know 


Qn- 
cat = age epee: 
Qn-2 


Since dy = 4n4n-1 + Qn-2 We find 


Qn Qn-2 1 
Qn-1 Qn-1 ee eee 


and our assertion follows by induction. 


1, §2. THE CONTINUED FRACTION OF A REAL NUMBER 


Consider first briefly the special case of a rational number a. Let ay be 
the largest integer < «. If a is not an integer, we can write 


1 
&=Apo +— 
hy 


with a, > 1, and a, is again rational. Inductively, we let 


Oy = a, AF > 
n+1 


where a, is the largest integer <a, and «,,, > 1. We can do this pro- 
vided «, is not an integer. However, our process will stop after a finite 
number of steps. Indeed, if «, = a/b is a rational number, not an integer, 
with positive integers a, b, then 


a—ba, c 
a, — a, = i 


b b 
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with c < b. Then 
b 


Cn+1 = 7 


and hence the denominator of «,,, is smaller than the denominator of a. 
So the process stops, and we can write our rational number « in the 
form 


g 


O= (do; ay, <-+-0, |, 


with integers a; (i =0,...,n), and a;=1 for i= 1. Observe that we have 
a choice for the last partial quotient a,, namely we can write « in the 
above form with a, equal to an integer > 1, or also in the form 


& = [do, a,,...,a, — 1, 1]. 


Thus the length of the continued fraction expansion of a rational number 
may be taken to be either even or odd. 

Let « be a real irrational number. We can determine a continued 
fraction for a by writing as before a =a) +1/a, with ay equal to the 
largest integer < a (that is, ag = [a]), and a, > 1. Inductively, we let 


a, = a, ak ’ 
On+1 


where a, is the largest integer <«,, and a,,, > 1. Since « is irrational, 
the sequence a,,%,,... does not terminate. We have, in the notation 
of §1, 

a= [ao, Ay, +++5An, tn+1 J 


for any n= 0. We shall also write symbolically the infinite expression 


Coa 


From Corollary 2 of Theorem 2, we obtain a sequence of relatively prime 
integers p,,4, With g, = 1, belonging to the continued fraction [do,...,a,], 
and thus the relation 


Pn 
a [ao, words 


is now a relation between real numbers, not any more between indeter- 
minates. Furthermore, p,/q, is a reduced fraction, which will be called 
the n-th principal convergent of «. We call a, the n-th partial quotient 
of a. 
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The formalism of §1 now applies to the continued fraction for a. For 
instance, the Corollaries 3 of Theorems 2 and 3 must now be written 


(- iy 
n+ 24n+1 a 3 Qn 


Gn+1% — Pati = 


and 


(—1)"a,,42 


.|), = 7 
a On+29n+1 2 Qn 


We always have 
a, <a, <a, + 1, 


and a, 21 for all n=1. Hence the denominators q, are all positive 
integers, and form an increasing sequence, 


O< 41 <°"° <4, < Gat <°**- 


Theorem 5. For even n, the n-th principal convergents of « form a 
strictly increasing sequence converging to a. For odd n, the n-th princi- 
pal convergents of « form a strictly decreasing sequence converging to a. 
Furthermore, we have 


1 1 
< <|q,% — Dal < : 
Ps hei, RO rer? | Qn+1 


Proof. The first assertion follows from Corollary 2 of Theorem 3, §1, 


and Corollary 1 of Theorem 2. So does the inequality on the right. The 
left inequality follows from Corollary 1 of Theorem 3, §1, namely 


Pn+2 Pn ike. 


Qn+2 Dn 


Gn+2 Gn+2 


GWn+2On (Qn4249n+1 “Tr An) In 


We divide numerator and denominator by a,,, and use the fact that 
Qn+2 2 1 to conclude the proof. 


The picture illustrating Theorem 5 may be drawn as follows: 


pt tt 


P2m-2 P2m eae ee P2m+1 Pom-1 


d2m—2 dam Gam+1 G2m-1 
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Since 4,4; >4,, we find that our convergents give us solutions of the 
inequality (1), namely 


1 
land = Pr Ss 
q 


We shall determine in §3 what other possible solutions may exist. 
We observe that for n = 1 we have 


1Gn 3 Pal = nll. 
Corollary. For n = 2 we have 


dn—1%ll = QnllQn@ll + Wdn+1 ll, 
whence 


anal < Ildn—1 ll, 


_ Maes 
a, =| Men 


Proof. We use Theorem 1 to express g,4;% — Pn+,, and then use the 
fact that q,4:% —P,+, and q,%—p, have opposite signs by Theorem 5. 
This proves the first relation of the corollary. The others are immediate 
consequences. | 


and 


We shail now characterize the principal convergents to « by an order- 
ing property. 
A best approximation to « is a fraction p/q (q > 0) such that 


llqall =|qx—pl, and = |iq’al| > Iga 


if 1 < q’ <q. Observe that the fraction p/q is necessarily reduced (i.e. p, q 
must be relatively prime) if it is a best approximation to a, for otherwise, 
we can write p = p’r, q = q'r with r > 1, and q’ <q, so that 


|g’ — p’| <|qa — pl, 
which is impossible. 
Theorem 6. The best approximations to « are the principal convergents 
to a. In fact, for n=1, q, is the smallest integer q>q,-, such that 
I]qoll < lldn—1 all. 
Proof. Let us first show that a best approximation is a convergent. 
Let a/b be a reduced fraction, b > 0, which is a best approximation to «. 


We must show that a/b = p,/q, for some n. Suppose that a/b < po/qo = ao. 
Then 


a 
la — ap| < a—$| <\ba— al 
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contradicting the hypothesis. Suppose that a/b > p,/q,. Then 


a p 
>|5—m 


2 
b = 


1 
Bi 


whence 


{ ii 
|ba — al > — =— 2 |a — aol, 
Qi Ay 


again contradicting the hypothesis. Finally, suppose that a/b lies between 
Pn—1/Qn-1 290 Pyi1/dn4,, but is not equal to either of these fractions. 
Then 

1 


GnQn-1 


a _ Pn-1 


Pr Pn-1 
< — = 
b Qn-1 


Qn Qn-1 


1 
aS 
bqn-1 


Hence q, <b. On the other hand, 


Ngee Pete 202 | ee 
bGn+1 Qn+1 b b 
whence 
lan% — Pal < = |ba — al, 
nt+i 


contradiction. This proves the first half of the theorem. 

We shall prove the converse, by induction on n. First for n = 0, since 
qo = 1, there is no q such that 1 < q < qg. Hence the definition of best 
approximation is vacuously satisfied by po/qg. Assume now that our 
assertion has been proved for p,/q, with n 20. We wish to prove that 
Pn+i1/Qn+1 18 a best approximation. Let q be the smallest integer > q, 
such that 


llall < lana, 


and let p be such that |/qa|| =|qa—p|. Then by the inductive prop- 
erty that p,/q, is a best approximation, we conclude that p/q is a best 
approximation also, and hence must be a principal convergent by what 
has already been shown. Since q is chosen smallest > gq, such that 
|qa\| < ||q,%l|, it follows that q = q,4,. But then p = p,,, (trivially), there- 
by proving our theorem. 


Corollary 1. If p/q is a principal convergent to «, and m is an integer 
with 1 <m < q, then 


1! 
at || ma. 
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Proof. Suppose that q = q,. By Theorem 5, we have 


1 
ao = 
9g, < Nn 


and by Theorem 6, ||q,_,«|| < ||ma||, as was to be shown. 
Corollary 2. If a/b is a reduced fraction, b > 0, such that 


2 
oe 


a 
04 —_ — 

b 
then a/b is a principal convergent to «. 


Proof. By the theorem, it will suffice to prove that a/b is a best 
approximation to «. Let c/d be a fraction, d > 0, c/d # a/b, such that 


1 
Pe cleans — al. 
|dx —c| < |ba al<—, 
Then 
1 c h(a : fe a 1 | b+d 
ile © 8) Sy = C2 ee ae Sa 
Ba bla |" gol Bl fea 2b? > oe 


From this we conclude b < d, whence a/b is a best approximation to a, 
as was to be shown. 


I, §3. EQUIVALENT NUMBERS 


a b 

Diced 
with integral components a, b, c,d having determinant +1 (ie. ad — be = 1 
or —1) is in fact a group, for the product of two such matrices and the 
inverse of such a matrix again have determinant +1. Let G be this 


group. If « is an irrational number, and o is as above an element of G, 
we define 


The set of matrices 


ax+b 
oa = : 
cat+d 


Then one verifies by brute force that if 0, t¢ G then 


a(ta) = (ot)a and la=a 
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if I denotes the unit 2 x 2 matrix. Thus G operates on the set of 
irrational numbers, and we shall say that two irrational numbers «, 8 are 
equivalent if there exists o€G such that ox =. It is trivially verified 
that this is an equivalence relation. 
Example 1. For « irrational, we can write 
Wide aso, Geno: 


By Theorem 1 of §1, and its corollary, we obtain for n 2 1, 


c= Pn-1%n a Pn-2 
Qn-1 On a n-2 


a ee fe 
Ci = 
Qn-1 Fn-2 


and call o,_, the (n — 1)-th continued transformation of «. We see that 


Let 


4= O,,-1 oA 
Furthermore, Theorem 2 of §1 shows that o,_, is an element of our 


group G. Thus « is equivalent to «, for n21, and consequently all 
numbers «, (n = 1, 2,...) are equivalent to each other. We also note that 


if we let 
aw. Il 
A,={ ” 
BG 


then det A, = —1, and from Theorem 1 of §1, we find 


In the next theorem, we give a characterization of the situation 
described in our example. 


Theorem 7. Let «, B be irrational, and 


_{a b ap + b 
o=(° eG & = OB — ee: 


Assume that B > 1, and c>d>0. Then b/d and a/c are two successive 
principal convergents to a, Say Py—-2/qn—2 and P,-1/q,-1, and B = d,. 
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Proof. Note that a,c are relatively prime because ad — bc = +1. We 
can express a/c as a continued fraction, 


a ; Ph-i 
~ = [do,...,4,-1] = ——, 
c : Qn-1 


and we have a = p,_,, C= 4q,-,;._ We may choose n so that 


Ph-19n-2 =e Up Ue) = Ss: 


where € = ad — bc (because the continued fraction of a rational number 
can be shortened or lengthened by 1 artificially). Since 


ad — bc = p,-,d — q,-,b =€ 
we find 


Daail@eG, >) = Gy-a() — p,-5). 


Since p,-1, ,-; are relatively prime, it follows that g,_, divides (d — q,-_>). 
But q,-2 <= 4,-; and d <q,_,. Hence 


|2 — q,-21 < Gn-15 
and therefore d — q,_, = 0. Then b — p,_, =0. Thus we can write 


3 Pe PP aes 


_ Gn-1B 2m Gs. 
This means that 


e= [ao, eee sAn—-15 Pl. 


Since B > 1, it follows that the above expression is the continued fraction 
expansion of a, and that B = «,. We then see finally that a/c and b/d are 
consecutive principal convergents, as was to be shown. 


Theorem 8 (Serret). Let «, 8 be irrational numbers. They are equivalent 
if and only if «,=f,, for some pair of integers n, m2 1, or equiva- 


lently, in their continued fractions 


a= [a, a,,42, dul: 
B — [do; ie b,, aly 
we have a, = b,., for some | and all n sufficiently large. 


Proof. Assume that there exist integers k, / 2 1 such that «, = f;, Le. 


“g= [ao, ay; vos, Ay_35 a, |, 


B =H [bo, D,, eee 1 Bl, 
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and a, = f,. Since we have seen that « is equivalent to %,, and f is 
equivalent to f,, it follows that « is equivalent to B. 
Conversely, assume that «, 8 are equivalent, say 


eee 
~ cat+d 


B 


(eyes 


with ad—be=+1. Without loss of generality, we may assume that 
ca + d > 0 (otherwise, replace a, b, c,d by their negatives). Let o,_, be as 
in Example 1, so that a =o,_,¢,. Then 


B = 00,,-1%> 


a) ee 
etre CP j24 + ddyai CPp=aet 44,2 Ca 


CPn-1 ae dqy—1 = Qn-1 (< a 2 ‘) c’, 


n-1 


and 


We have: 


CPn-2 + ddy-2 = Qn-2 (< ae 4° ‘) =’. 
n-2 


We take n large, so that p,_,/q,-; and p,—2/q,-2 are close to «. Then c’ 
and d’ are both > 0, and also a, > 1. Finally, we can take the parity of 
nso that c'’ >d’. Then all the conditions of Theorem 7 are satisfied, and 
we conclude that «, = B,, for some m. This proves our theorem. 


Examples. Suppose first that 


a= lage PO Passed 


Then it is easily verified that 


~~ | eee if'a, >. 
[—a — 1, a, + 1, a3, a4,...] ed) 


It is also easy to take the inverse. We have from the definitions: 


oo if «> 1, 
[atvas assed i Olecore Ie 
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I, §4. INTERMEDIATE CONVERGENTS 


We continue with more results of Lagrange, and define fractions lying 
between the principal convergents of an irrational number a. We first 
note that if a, b, c,d are non-zero real numbers such that 


The proof is trivial by cross multiplying. Furthermore, if r < s then 


(@h S> TRE “ar Se 


meer pasa’ 


in other words, the quotient (a + rc)/(b + rd) is a strictly increasing func- 
tion of r. For any integer r = 0 we let 


Par = YPn+1 a2 Pn and Qn,r a TQn+1 ae Qn: 


If r=0, then p,o = p,, and if r=a,,. then p,,=DPa+2- Similarly for 
Gn,r- We shall be interested in the values of r such that 


1 S i An+2> 
and call the fractions 


Pn,r = TPn+1 a Pn 
Qn,r l'On+1 aS a, 


S70. 1; 


the n-th intermediate convergents of the continued fraction [dp, 4,,...], 
or of a if this continued fraction is the one associated with «. An inter- 
mediate convergent or a principal convergent will be called a convergent. 

We note that the denominators of the intermediate convergents and 
the convergents form a strictly increasing sequence 


i“ S Ces ae S Clap S Cys Se < n+2 ne 


16 GENERAL FORMALISM [T, §4] 


Theorem 9. For n even, we have a strictly increasing sequence 


+ 
ohn e.g mt. a 
Qn Qn.r  Anyr+1 Qn+2 


KG NE 


and a similar decreasing sequence for n odd. Furthermore, 


Qn,r+i1Pn,r —~ Pn,r+14n,r = ee ‘eg ° 

Proof. The increasing sequence comes from the remarks on inequalities 
made at the beginning of the section. The last relation follows from a 
trivial computation, using Theorem 2, §1. 

We conclude that p,, and q,, are relatively prime, and thus the 
intermediate convergent p, ,/q,,, is in reduced form. 

The next result, which is important for the determination of the solu- 
tions of the fundamental inequality |qa — p| < 1/q is due to Grace [10] 
(cf. also Adams [1]). 


Theorem 10. If p,q are non-zero integers, q > 0, satisfying the inequal- 
ity |« — p/q| < 1/q?, then p/q is a convergent of «, and is in fact equal 
to some Dy.p/Qn,, With r = 0, or r = 1, or r=a,42 — 1. 


Proof. For the first statement, assume, say that « < p/q, the other case 
being proved in a similar manner. If p/q is not a convergent, then there 
exist two successive convergents P/Q and P’/Q’, such that 


a < P/Q <p/q< P/Q 
and P’Q — PQ’ = 1, by Theorem 9. Thus 


i 
2) 


a 
qv 4q q 


J 
qQ 
ee ae ae 
Oq QO q @ Q QQ 


and 


These estimates are contradictory, and our first assertion is proved. 


Conversely, we must determine which convergents Pn.r/Qn,r Satisfy the 
desired inequality. We need a lemma. 


Lemma. If p,,,/q,,, is an intermediate convergent of a, then 


_ ("G2 — 1) 


Qn,r% — Dy, : 
ie ee Oar Qn 
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Proof. We have the continued fraction 


& = [o, G1, ..+,4n415 On+21- 


The relation of Corollary 3 of Theorem 3, §1, is a relation between 
indeterminates. Consequently it holds for the special values of the pre- 
ceding continued fraction, and we find 


(= 1)"on+2 


Gn% — Py = ———-_.. 
alee On+29n+1 + An 


Similarly, using Corollary 3 of Theorem 2, §1, we find 


(— Lys 


Qnt1% — Pati = —————: 
oa On+29n4+1 ale An 


Multiplying this last expression by r and adding the preceding one, we 
find the equality stated in the assertion of the lemma. 


We observe that r < «4, because r<a,,,— 1, and consequently that 


G42 —? 


da. — oe) = 
a : On+29n+1 at qn 


We must find a necessary condition on r for this expression to be 
< 1/q,,,, in other words, 


Xn+2 = r PE 1 
Ont+29n+1 + 4n "On+1 + Qn 


A trivial manipulation shows that this inequality is equivalent to 


Suppose r > 0. Since r < «,42, this inequality implies 


r(1 Payor )< 1M 
Xn+2 


Suppose r < a,,, —1 and hence r <«,,2 — 2 because r is an integer and 
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42 is not. Then 


whence r < a@,,,/2. Consequently 


led 
ae sw.) a 
( 2%n+2 


and r < 2, so that r= 1. This proves Theorem 10. 


HI, §4] 


We can generalize to intermediate convergents one of the inequalities 


proved previously for the principal convergents of «. 


Theorem 11. Let p,,,/q,,, be a convergent of « with 


0 = Ti < do a 1. 
Then 
1 


nr ag Qn,r+1 


Proof. By the lemma of Theorem 10, we must check if 


1 . Cnies —T 
(2r tr N4n+1 a 24n On+29n+1 + In 


This inequality is equivalent with 


1<(1- Var 4 ja 
OCn+4 


and is implied by 


1s(1- "Jar+ 0) 


On+1 


<< Gare a Pix: 


The same kind of trick used in the proof of Theorem 10 now shows that 


this inequality is satisfied for 0 <r < a,,, — 1, as was to be shown. 


A solution p/q with relatively prime integers p and q, q=1, of the 
inequality |q« — p| < w(q)/q, with some positive function «, will be called 
an w-convergent of aw. The 1-convergents of « are therefore the reduced 
fractions p/q, with q = 1, of the fundamental inequality |qa — p| < 1/q. 
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We order the 1-convergents by increasing denominators. By Theorem 
10, the sequence of denominators looks like 


aS Qn+1 < Qn,1(?) x cae (?) < Qn+2 ss 
and the question marks mean that the intermediate terms q,,, or Ga 
may or may not be present. 


As a special case of Theorem 11, we find: 


Corollary. If p/q, and p'/q’ are two successive 1-convergents of a, then 


7 = ga — pi. 
q+q 


This section more or less concludes our study of the formalism of 
continued fractions. For a discussion of more specialized topics, cf. 
Perron’s excellent book [21]. It is a problem to extend the results of 
this chapter to simultaneous approximations, those being described by 
investigating 

qi X4 a eke GaAs 


where X; is a vector in some higher dimensional space, and the norm is 
that of the torus. A first attempt was made by Perron [22], and was 
recently pursued by Bernstein [5], [6], [7]. For a description of possible 
eventual applications to contexts of algebraic geometry, cf. [17], [18]. 


CHAPTER Il 


Asymptotic Approximations 


ll, §1. DISTRIBUTION OF THE CONVERGENTS 


We begin by an old result of Dirichlet. 


Theorem 1. Let « be a real number, and N a positive integer. There 
exists an integer q,0<q<N such that |\qa|| < 1/N. 


Proof. Cut up the interval [0, 1] into N equal segments of length 1/N, 
and consider the N + 1 numbers 


Ow, la, 20, ...,Na 


modulo Z. Two of them must lie in the same segment (mod Z), say ra 
and sa with r<s. We let g=s—r, and obtain 


1e4) 
laall< > S-. 
q 


as desired. 


We are interested in a lower bound for the integer g of Theorem 1. 
Let a be an irrational number. Let g be a positive function, which will 
always be assumed to be increasing (not necessarily strictly), and > 1. 
We shall say that « is of type <q if for all sufficiently large numbers B, 
there exists a solution in relatively prime integers q, p of the inequalities 


laa—p|<1/q and B/g(B)<Sq<B. 
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Theorem 2. Let {p,/q,} be the sequence of principal convergents to a, 


and let f be an increasing function 21 such that for all n sufficiently 
large, 


1 
—.—~ S$ |q,% — p,|. 
ay 
Then « is of type < f. 


Proof. Since |q,% — p,| < 1/9,4;, we conclude that 


Gee <n Gs). 


Given N large, we find n such that q, < N Sq,i,. Then 


N .N 
eee < gy ON 


F(N) ~ £(4n) 
thereby proving that « is of type < f. 


Theorem 2 admits a partial converse, which shows that a type for a 
determines some kind of lower bound for |gqa — p| with q, p relatively 
prime. 


Theorem 3. Let « be of type <g. Assume that the function t/g(t) is 
strictly increasing, and let g* be its inverse function. Then for any 
sufficiently large integral solution q, p of |qa — p| < 1/q, with q, p rela- 
tively prime, we have 
1 
—eur See 
q+ 9*(q) 
Proof. Let p/q be a 1-convergent of a, and let p’/q’ be the 1-conver- 

gent of « with smallest denominator > q. Then by hypothesis, 


, 


q 
g(q') 


Sq<4q 
whence q’ < g*(q). By the Corollary of Theorem 11, Chapter I, §4, we 
conclude that 

1 


= = go — pi, 
Go \¢) aed 4.4 


as was to be shown. 


Remark 1. Suppose that g is a function which grows reasonably 
slowly, namely such that there exists a constant c 2 1 for which 


g*(t) S ctg(t) 


22 ASYMPTOTIC APPROXIMATIONS : [II, §1] 


for all t sufficiently large. Then we can rewrite our inequality 


ee 
ia P| 


This happens when g is constant, or grows like the log, or a power of 
the log. When g(t) = t‘, then one gets a slightly different estimate. 


Remark 2. In Theorem 3, and also in other applications (viz. Chapter 
III, §2) it is more useful to deal with a variation of the notion of type. 
Thus we may say that « is of cotype <g if given a sufficiently large 
number B, there exists a solution in relatively prime integers q, p of the 
inequalities 

lqx—pl<l/q and B<qs Bg(B). 


As we saw in Remark 1, in most applications a type and cotype can be 
taken as the same function. 


To get some idea of possible types for numbers, we shall now prove a 
simple theorem of Khintchine. We recall that a set of numbers is said to 
have measure 0 if given « > 0, the set can be covered by a countable 
number of intervals, such that the sum of the lengths of these intervals 
is <€. 


Theorem 4. Let w be a positive function such that 


Ms 


w(q) 


q=1 


converges. Then for almost all numbers « (i.e. outside a set of measure 
0), there is only a finite number of solutions to the inequality 


lal] < W(q). 


Proof. Given € > 0, select gg such that 


y W(q) <€«. 


9240 


We may restrict our attention to those numbers « lying in the interval 
[0, 1]. Consider those for which the inequality has infinitely many solu- 
tions. For each q 2 qo, consider the intervals of radius w(q)/q surrounding 
the rational numbers 
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Every one of our « will lie in one of these intervals because for such a 

we have 

eae 
q 


q 


The measure of the union of these intervals is bounded by the sum 


yy ge <2 


4540 


as was to be shown. 

For example, we can take w(q) = 1/q(log q)'** for any e>0. Thus 
we can take the function f(t) = (logt)'** for almost all numbers in 
Theorem 2. 


Theorem 5. Let w be a positive function such that 
» wa) 
q=1 


diverges. Then for almost all numbers a, there exist infinitely many 
solutions to the inequality ||qa|| < w(q). 


We refer to Khintchine’s book for the proof of Theorem 5. 
Theorems 4 and 5 will be called Khintchine’s convergence and divergence 
theorems respectively. 


li, §2. NUMBERS OF CONSTANT TYPE 


There is a special kind of numbers which provides useful examples, and 
is especially easy to work with. They are characterized by the properties 
of the next theorem. 


Theorem 6. The following properties concerning an irrational number a 
are equivalent. 


CT 1. There exists a constant c >0 such that for all integers q >0 
we have ||qa|| > c/q. 

CT 2. For any positive function W with convergent sum > W(q), the 
inequality 


Ilqa\| < W(q) 


has only a finite number of solutions. 
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CT 3. There exists a constant c>0 such that, given a sufficiently 
large integer N, there exists a relatively prime solution q, p 
of the inequality |qa — p|<1/q, and N<q<cN. (In other 
words, « is of constant type.) 

CT 4. If [ao,@,,42,...] is the continued fraction of a, then there 
exists a constant c > 0 such that a, <c for all n. 


Proof. Assume CT 1, and suppose that yw is a function such that 
|qa|| < W(q) has infinitely many solutions. Then c/q < W(q) for infinitely 
many gq. We contend that the sum ) y(q) diverges. 

Dividing w by c, we may assume without loss of generality that c = 1. 
Let q, <q, <°': be the increasing sequence of q such that W(q,) > 1/q,. 
Define (q) = 1/q, for q,-1 <4 <4,-. Then go <¥y, and it suffices to prove 
that }° p(q) diverges. Then 


1 1 
Yu@2=> e(q@z a + (q2 - 41) —+ (43 — 42) — + 
qQ1 qQ2 q3 


Take n=n, large. The first n terms of this series have a lower bound 
given by 


1 14g 
aaa + (Gn — dns) = 


Thus for n large, we get a contribution > 4 to our sum. We repeat this 
procedure with a number n, which will give a contribution greater than 


to our sum, and so on with n;,.... In this manner, we see that the sum 
diverges, and CT 2 is proved. 

Assume CT 2. We shall prove that « satisfies CT 1 by an argument 
due to Schanuel. Suppose that « does not satisfy CT 1. Then we can 
find a sequence of integers q; with 


1 <qi9g3 <= 


such that ||q;«|| < 1/2'g;. Let 


e7t/4i, 


W(t) = 


2 


ée 
1 2'q; 


Then w(q;) > 1/2!q; for j = 1,2,... and the sum for w converges. This is 
a contradiction, which proves that « satisfies CT 1. 
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We observe that Schanuel’s function is very smooth, and behaves as 
well as possible from the point of view of convexity. Thus if CT 2 is 
assumed only for such functions, it still follows that a satisfies CT 1. 

The equivalence of CT 1 and CT 3 is a special case of Theorems 2 and 
3. The equivalence of these with CT 4 follows from the fact that at most 
two n-th intermediate convergents are also 1-convergents, by Theorem 7 
of Chapter I, §3. This proves our theorem. 


Numbers of constant type are also said to have bounded partial quo- 
tients, in view of CT 4. 


Example. Let D be a positive integer which is not a square, and let 
a=at+ b./D where a, b are integers. Then a is of constant type. This is 
trivially seen as follows. Suppose that |qa — p| is small, so that « — p/q is 
small. Let «’ =a— b./D be the conjugate of «. Since p/q approximates 
a very closely, we conclude that a’ — p/q is approximately equal to a’ — a. 
But (qa — p)(qa’ — p) is a non-zero integer, of absolute value = 1. If 
|qa — p| < c/q for some small c > 0, then 


|qax’ — p| 2 q/c. 


However, qa’ — p is approximately equal to g(a’ — «). This shows that c 
cannot be arbitrarily small. 


Example. If « is of constant type, and m/n is a rational number # 0, 
then ma/n is also of constant type. The easy proof will be left as an 
exercise to the reader. 


In view of Khintchine’s divergence theorem, we see that given an 
integer n>0, the set of numbers a for which there is only a finite 
number of solutions of the inequality ||qa|| < 1/nq has measure 0. Call 
this set S,. If m>n then S,cS,,. Every element of S, is of constant 
type, and conversely, every number of constant type lies in some S,. 
Since the countable union of sets of measure 0 also has measure 0, it 
follows that the numbers of constant type form a set of measure 0. 

No simple example of numbers of constant type, other than the one 
given above, is known. The best guess is that there are no other 
“natural” examples. 
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Throughout this section, we let w be a positive function <1, decreasing, 
such that 


2 v(@) 
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diverges. We let 7 
VN} = | y(t) dt. 
it 


For each positive integer N and irrational number «, let 27 ,(N) be the 
number of solutions in integers q, p of the inequalities 


0< qa—p<w(q) and 1 SN. 


To simplify the notation, we shall omit the + sign, and also we usually 
omit the indices «, y on 4. It is natural to ask for an asymptotic esti- 
mate for A, but such an estimate was proved only recently for almost all 
numbers. We shall state this result (without proof). We first recall some 
terminology. 

If F,G are two functions of a real variable, and G is positive, we say 
that they are asymptotic and write F ~ G if 


lim F(x)/G(x) = 1. 


We say that F =O(G) if there exists a constant C>0O such that 
|F(x)| < CG(x) for all x sufficiently large. We say that F = o(G) if 


lim F(x)/G(x) = 0. 


xo 


Theorem 7. For almost all numbers «a, we have 
A(N) = P(N) + o( ¥(N)). 


A special case of Theorem 7 was first stated by Leveque [19]. The 
general theorem was proved by Erdés [8] and Schmidt [24]. In this 
book, we are principally interested in specific numbers, and we shall omit 
the proof of Theorem 7, but give a partial result (Corollary 3 of Theorem 
8 below) consistent with our point of view. We point out, however, that 
Schmidt obtains further important generalizations, e.g. higher dimensional 
ones, and also has a very good error term. This is important, because in 
dealing with specific numbers, the expression of the error term reflects 
the special nature of the number under consideration in an essential way. 
For further work on this, cf. also Gallagher [9]. 

It is a problem to determine specific numbers, and functions w for 
which 4 has a similar asymptotic property. For the statement of the next 
theorem, we introduce some notation. We write f >g and say that f is 
much larger than g if there exists a positive function h tending to infinity 
such that f = gh. We also say that g is much smaller than /,. 
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Theorem 8. Let a be an irrational number of type <g. Write 
W(t) = w(t)/t. Assume that w>g, that w is increasing to infinity, and 
that w(t)"?g(t)'?/t is decreasing for all t sufficiently large. Then 


A(N) = P(N) + of [seh eo" a) 


1 


Remark. If 7 is a function tending to 0, then one verifies easily that 


for N > 00, 
| “on | 


ry t = o('P(N)). 


Consequently, since @ >g, we see that the error term given in the theo- 
rem implies the asymptotic result 4 ~ YP. 


If w is a number such that |qa — p| > 1/qf(q) for some increasing 
function f and q, p relatively prime, then we know by Theorem 2 of §1 
that we can take g =f, whence the asymptotic result holds whenever 
a> f. We. have two interesting special cases: 


Corollary 1. If « is of constant type, then 


a(N) = P(N) + o( | oO at] 


for any function @ > 1. 


Corollary 2. Let 0<a<1 and let w(t) = at. Then A(N) is the number 
of pairs of integers q, p satisfying 


O0<qua—p<a and l<q<N. 


aN) = an + of |" 0h ae). 


1 


We have 


The error term is o(N) if g(t)/t tends to 0 ast. 


When @ is as in Corollary 2, then the problem of estimating / is 
known as the equidistribution problem. It determines the number of 
integers q such that qa (mod Z) lies in the interval [0,a], satisfying 
1<q<N. When A(N) is asymptotic to aN, we interpret this as saying 
that the numbers qa (mod Z) are equidistributed. Corollary 2 determines 
the connection between this equidistribution problem and the type of the 
number a, by means of the error term. This particular case had been 
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considered long ago, notably by Weyl [29], and in a manner more 
closely related to the point of view taken here, by Ostrowski [20], and 
Behnke [4]. Instead of working with the type as we have defined it, 
however, these last-mentioned authors worked with a less efficient way of 
determining the approximation behavior of « with respect to p/q, whence 
followed weaker results and more complicated proofs. The function a, 
which shows itself to be quite important in the present estimates was 
introduced in [15]. 

Theorem 8 also implies a statement about almost all numbers, since 
we can apply Theorem 4, §1, to these. If g(t) =(logt)'** then the 
Khintchine convergence theorem implies that almost all numbers are of 
type = g. Thus: 


Corollary 3. Let w be a positive function such that w >log'**. Then 
for almost all numbers « (the exceptions being on a set of measure 0, 
depending on w), we have 1, , ~ Y. 


The proof of Theorem 8 will involve first a special case of Corollary 2, 
as in Lemma 1 below. 

It is convenient to abbreviate the remainder of a number € (mod Z) 
between 0 and 1 by R(é). Thus R(é) is the unique number € — p (with 
some integer p) such that 0 < €—p< 1. 


Lemma 1. Let 0<a<1. Let p,q be relatively prime integers, such 
that |qx—p|<1/q. The number of integers n among q consecutive 
integers such that R(na) <a is equal to ga + O(1). 


Proof. Write 
a : with 6] <1. 


Let n =n, + v, with v= 1,...,...,q. Then 


m) 
no = (Mo + va = now + va = now +P + o. 
q 4 


The rational numbers vp/q (v = 1,...,...,q) are equal to 


q 


0 —1 
—,-,...,——— (mod Z), 
q q 


1 
q. 


up to a permutation. The error vd/q? is bounded by 1/g. We can write 
Noa =r/q+e for some integer r, O<r<q and |e|<1/g. Hence the 
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numbers na (mod Z) are precisely the numbers 


r+yv 
q 


+e, (mod Z) 


with |e,| < 2/q. Up to a permutation, these numbers are nothing but the 
numbers 


: +é€, (mod Z) 


with 0 < » <q and |e,| < 2/g. Now we have 


B L 
R Hee jatre 
G ae 


except possibly when p = q—2, or » <2 which occurs for at most five 
values of yu. Thus the number of desired integers n is equal to the 
number of integers u with 0 < p < q such that 


m 
—~+€,<a, 
q H 


up to a bounded error term. The number of solutions of this inequality 
is bounded from above by the number of solutions of 


2 
eae 
qq 


or equivalently, u < qa+ 2, which differs from qa by a bounded error 
term. Similarly, we obtain a lower bound which also differs from qa by a 
bounded error term, and thereby prove our lemma. 


Lemma 2. For all N sufficiently large (depending on w) and all integers 
q, p relatively prime, q > 0, satisfying the inequalities 


Ng(N)*? 
= < pia dh alae 
0 < |qa — p| < I/q and lsq< oN)” 
we have 
; N N 1/2 N 1/2 
AN) - 20 -9)= | wie) de + 04% y+ 41, 
N-q 


where |@| < 4, |@,| <c,, and c, is an absolute constant. 
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Proof. We note that 4(N) — A(N —q) is the number of integers n 
satisfying 


0 < R(na) < W(n) and N—qsn<vN. 
Let 
_ a('"g(0"” 
— 


E(t) 
We contend that for N sufficiently large, 
0<¥(N — gq) — WN) S 2E(N). 


To see this, note that 


o(N — 4) oN) 
N-@q N 


0S ¥(N — q)— ¥(N) = 


_ No(N — 4) — (N — go(N) 
N(N — q) 


We replace w(N —q) by w(N), making the right-hand side bigger. 
Similarly, we can replace N—q in the denominator by N/2, because 
g(t)/w(t) +0 as too. Using the assumption that q < Ng(N)*?/a@(N)*?, 
we see that our contention follows at once. 

We conclude that 


W(N) S W(n) S WN — gq) < WN) + 2E(N). 


We can determine bounds for A(N)—A(N —4q) replacing w(n) by 
W(N) + E(N) in the inequality 0 < R(na) < w(n). By Lemma 1, we obtain 


A(N) — AN — q) = qu(N) + 2gE(N) + O(1). 


On the other hand, since y is decreasing, 


N 


qu(N) Ss le W(t) dt < qu(N — q) S qu(N) + 2gE(N). 


Hence 


N 
A(N) — A(N — q) = | y(t) dt + OqE(N) + 6,, 
N-q 


thereby proving Lemma 2. 
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We may now give the main part of the proof. For N sufficiently large 
let 


a 


p= Nainy? 
a w(N)*2 = ’ 


and select qg,p such that 0 <|qa —p|<1/q and B/g(B)<q<B. Then 
trivially, since g(B) < g(N), 


Ng(N)"? 
@(N)"? 


N 
Tea a 
o(N)g(Ny? = 4 


and 
0 < |qa — p| < 1/q. 


Using the fact that o(t)'g(t)"?/t is decreasing, and the left inequality for 
q, we get 


Eel. 


{ (790? ay aenirat 


N-q t _ 


By Lemma :2, it follows that 


A(N) — A(N — q) = [ Wd) dt + Oy, lo o( ay” a 


N-q t 


with |6y ,| Sc, +5. Repeating our argument with N — q instead of N, 
and taking the sum inductively, we find that 


A(N) = YN) +6 | 5 OE 5. O(1), 


1 


with |6| <c, +5. This proves our theorem. 


Aside from possible quantitative generalizations to higher dimensions, 
one also faces the problem of describing the asymptotic behavior to the 
single number « when g grows slower than qm. It is then not necessarily 
true that A is asymptotic to (cf. Chapter IV and especially Chapter V). 
In the next section, we shall describe a method which may sometimes be 
applied to such cases. 

As for the higher dimensional case, one may begin by considering 
linear combinations 


eae = 4+ Geenll 


from the present point of view. Khintchine’s transference principle gives a 
weak relation between the size of the above combination, and simultane- 
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ous approximations to the numbers a,,...,...,4m, (which one assumes 
linearly independent over the rationals). (Cf. for instance, Cassell’s book 
Introduction to Diophantine Approximations.) It would be interesting to 
formulate a quantitative asymptotic transference principle, depending on 
a generalized notion of type for linear combinations of several numbers. 


, §4. RELATION WITH CONTINUED FRACTIONS 


We shall describe here in a formal setting the method used by Adams to 
determine asymptotic approximations to e and other numbers [1], [2]. 
A special case will be carried out in Chapter V. 

Let 


a = [do,a,,...] 


be irrational as usual. From the relation 


Qn+1 = An+1 qn a Qn-1 
we find 


1 1 ) 
+ ear | 1 ee 
dn is Qn/Qn-1 ty ( An+14y as Pn 


with some p, > 0. By induction, it follows that 


ROA me 
= Vow 


Let A(n) = a,°*-a, and let P(n) be the product 


n 1 
P(r) = |] (: + 1] 


A(n) = q, S A(n)P(). 


so that 


Suppose that there exist strictly increasing functions A, and A* such that 
A,(n) S A(n) and A(n)P(n) S A*(n). 

Then we obtain 

(x) A,(n) Sq, S A*(n). 


Let g, and g* be the inverse functions of A, and A*, respectively. 
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Let A)(N) be the number of solutions in relatively prime integers g, p 
of the inequalities 


(1) O<q«—p<ti1/q and OZ aan, 

and let A(N) be the number of solutions of these same inequalities with- 
out the restriction that q,p be relatively prime. We wish to find an 
expression for 4, and A in terms of the continued fraction for a We 


shall obtain it under the following: 


Assumption. The only relatively prime solutions q, p of the inequality 


qa — | a 
q 
for q sufficiently large are given by the principal convergents to «. 
Given N, let n be such that 
GN =< Gin. 


Then from (x) we find n< g,(N) and g*(N) < g*(q,41:) S$ n+ 1, or more 
clearly, 
g°(N) — 1 Simseg, UN): 


A principal convergent p,/q, will satisfy g,« — p, > 0 if and only if v is 


even, by Theorem 5 of Chapter I, §2. Hence by definition, 1,(q,) = 
4n + O(1), and we find the bounds 


(2) 


To find similar bounds for 2, we have but to count non-relatively 
prime solutions of our inequality for each v, 1 <v <n, and sum these. 
From Corollary 3 of Theorem 3, Chapter I, §1 (with «,,, instead of a,,., 
of course), we see that a positive integer k satisfies 


1 
|kq,a - kp, | << ka, 
if and only if 
k? < ee pees 
Oy+2 qy 


If a4, 22, 1<v<n, then trivially kq,<4q,,,. Hence the number of 
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integers k such that the multiples kq,,kp, yield a solution of (1) with 
ka, < dn+1 18 equal to 


ae + O(1). 


Hence we obtain the intermediate estimate 


n 


AGnts) = Lay” + O(n). 


< 


1 
veven 


o 


We can then determine /(N) easily from the inequalities 


Mn) S ACN) S AQn+1), 


using the bounds for n in terms of g*(N) and g,(N), namely: 


g*(N)—2 


94(N 
a,” — O(9,(N)) S A(N) S - ay!” + O(gAQN)). 


(3) 


v=1 y= 
eve veven 


All the terms appearing in the bounds (2) and (3) are expressed en- 
tirely in terms of the continued fraction for «, as was our goal. In the 
applications, the two functions A, and A* can be chosen such that g, 
and g* are quite close to each other, and in this way give an asymptotic 
estimate for A(N). 

The estimate of (3) is slightly coarse, due to the presence of the terms 
O(g,(N)). This is sufficient for the applications to continued fractions 
like those considered in Chapter V. One could get a somewhat more 
exact expression by being more careful in counting the k’s, satisfying the 
inequality 


i 1 1 
jo ee <44,+—+ 
qy Oy+2 a, Ay+2 


At present, one does not yet have sufficiently many examples of con- 
tinued fractions associated with classical numbers to be able to give 
useful and significant axiomatizations for this counting. 


CHAPTER _IIl 


Estimates of Averaging Sums 


Hl, §1. THE SUM OF THE REMAINDERS 


Let again « be an irrational real number. In considering the values 
R(nq), it is natural to form the sum 


3 R(na) 
n=1 


and to estimate its order of magnitude. Since one expects the values 
R(na) to be somewhat evenly distributed around 3, it is then better to 
investigate immediately the sum 


Sy = ¥: (Rina) — 3), 


which gives the average discrepancy between 5N and the sum of the 
remainders. We shall estimate this discrepancy in terms of a type for a. 


Theorem 1. Let « be of type < f, and assume that the function f(t)/t is 
decreasing. Then 

Dt 

1 
Proof. Let p,q be relatively prime integers such that 


1 
= and <q<N. 


ein 
f(N) 
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Then 
é 


; with {6| <1. 


aah 4 
q 


& 


We estimate the sum 


Sy — Sy-¢ = oe (R(na) — 4) 


see) 
= R| Na — v—-—v—]—<]. 
2 ( q gq) 2 
The fractions vp/q range over all fractions 0/q, 1/q,....(q — 1)/q and hence 
Sy — Sy-q = O(1). In particular, there exists an absolute constant c, such 


that 
N 
t 
5, - 5-<=0 | LO a 
N-q@ 


with |6|<c,, because the integral on the right satisfies the lower 


estimate 
[. te ele 
N-q 


Our theorem follows by induction, repeating our procedure replacing N 
by N — q. 


The preceding proof is that given in [18]. The sum of Theorem 1 had 
been considered classically, notably by Hardy—Littlewood [11], Ostrowski 
[20], Behnke [4], and Hecke [13]. The reader will find connections 
between this sum and other problems in these papers, especially as con- 
cerns the function defined by the Dirichlet series 


2? 


It would be quite interesting to investigate the analytic properties of this 
function of the complex variable s, and extend the results of Hardy- 
Littlewood and Hecke in this direction. This has to be done in connec- 
tion with special numbers, say algebraic numbers, or e, for instance. 

We observe that the error term of Theorem 1 is quite good. For 
example, if « is of constant type, then the error, i.e. Sy, is of the order of 
magnitude log N, which is quite small compared with the total number 
of terms in the sum. Similarly, we know that almost all numbers are of 
type < (log t)'**, so that for such numbers, the sum Sy is of the order of 
magnitude (log N)***, which is again quite small. 
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Ill, §2. THE SUM OF THE RECIPROCALS 


For subsequent applications, it is convenient to deal with a variation of 
the notion of type, and also to work only with the principal convergents 
to «. Let g be an increasing positive function = 1, and By a positive 
integer 210. We shall say that a is of principal cotype <g for all 
numbers 2 B, if given a number B > Bp, there exists a principal conver- 
gent p;/q; to « such that B <q, < Bg(B). In particular, if p/q and p’/q’ 
are two successive principal convergents to a, and q = Bo, then q’ < qg(q). 


Theorem 2. Let « be of principal cotype <g for all numbers = By. 
Then for all integers N = By we have 


N 1 
< 2N log N + 20Ng(N) + Ko, 
2 Rina) = og N + g(N) + Ko 
where 
Bog(Bo) 4 
Kes a=, 
°= 4A R(na) 


The same estimate holds if we replace R(na) by 1 — R(nq). 


Proof. The statement for 1 — R(na) follows just as for R(na). So we 
do just the first estimate. 


We shall need a lemma, estimating the sum of the reciprocals taken 
over certain consecutive integers. 


Lemma. Let p/q and p'/q' be two successive principal convergents to a. 
Let qo be an integer, 1S qo Sq. Let no be an integer 2 0, and assume 
No + do < qd. Then 


No+4do ‘ , 
Lil Rina) < q log q + 109’. 
n=No 


Proof. We write 


eta? isis, 
q 4 
and 
n=Nnot+v with v= l1,...,do. 
Then 


na = No + VA = Moa + ¥" +6, 


with |e,| < 1/q. The numbers ye (mod Z) are precisely the same as the 
¥ q 
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numbers 


F (mod 2), OSusq-1, 
q 
because p, q are relatively prime. Thus we can write 
a 1) 
noe = Not + 5 +e, (mod Z). 


(To each n we have associated a unique v and a unique p, which should 
be written v, and y,, but we omit the index n for simplicity.) We 
distinguish cases in estimating our sum. 

(a) n=q, or n is such that 


For such n, we see that R(nga + u/qg) occupies the position indicated by a 
cross in the following diagram: 


i ee eee q 
Qq qq @q q 


For such n, we use the fact that n < q’, and hence that 1/2q’ < R(na), by 
Corollary 1 of Theorem 6, Chapter I, §2, in other words, 


1 
ey 
R(na) ~ ae 


We observe that there are exactly 5 values of n in the present case, and 
the sum of the reciprocals for these n will be < 109’. 
(b) All other n, that is n 4 q, and 


Then R(na) = 2/q, and in this case the remainders 


R (nos t t) 
q 


lie in all the small intervals indicated by a cross in the following diagram: 
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Then 


R(noa + t) = R(noa = “i «) = R(na). 


We estimate 1/R(nx) replacing R(nx) by the fraction m/q lying immedi- 
ately to the left. If we take the sum of 1/R(na) for all n satisfying our 
condition (b), we see that this sum is 


q 
s¥ 7g 54108 4 


Adding the upper bounds in cases (a) and (b) yields an estimate which 
proves our lemma. 

We may now carry out the proper part of the proof of Theorem 2, by 
induction on N = By. Let p/q and p’/q' be successive principal conver- 
gents to a such that 

qeaN<q’. 


If q < Bo, then N < Bog(Bo) and we estimate our sum by Ko, trivially. 
Suppose q 2 By). Then 
q’ S$ 99(9) 


by definition of the cotype. We shall distinguish two cases. 


Case 1. N/2 <q. Then N < 2q, and we decompose our sum into two 
sums: 

N 
+>. 


qti 


le 


rMs 


We apply the lemma to each sum, and find 


1 
R(na 


using the monotonicity of log and g. 


< 2q log q + 20qg(q) S 2N log N + 20Ng(N), 


ote 


Case 2. q < N/2. From the lemma, we find 


N 


< 2g log q + 20 
weet Rina) = q log q qg(q). 


(Multiplying by 2 only makes the estimate worse.) By induction, 


— q) log(N — q) + 20(N — q)g(N — 4) + Ko. 


1 
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Using the monotonicity of g, we find 
(N — q)g(N — q) + 49(q) S$ Ng(N — 9) — ag(N — 9) + 9914) 
<Ng(N-—q) (because N—-q2q) 
< Ng(N), 


and similarly with log replacing g. This proves our theorem. 


Remark 1. The estimate of Theorem 2 is essentially best possible. 
Taking N = q—1 or N =q shows that both terms N log N and Ng(N) 
are necessary in the estimate. Furthermore, for « of constant type, one 
sees easily that the sum is 2 cN log N for some constant c > 0 and all N 
sufficiently large. 


Remark 2. An estimate for Ko itself can be given only in terms of g. 
We use the fact that 
N <q’ S Ng(N), 


and approximate a by p’/q’, rather than by p/q. Then we obtain the 
alternative estimate 


Y ss S 10Ng(N) log(Ng(N)) 


for all N 2 Bo. Thus sums in this chapter can be estimated only in 
terms of the type-—nothing else needs to be known about the number «. 


Remark 3. As usual, for almost all numbers we can take 
g(t) = (log t)'**, 


N61 
1 R(na) 


so that we obtain 


= O(N (log N)'**) 


for almost all numbers. Of course, the constant in O depends on ¢€ 
and «. 


Theorem 2 has various applications. We shall first use it to estimate 
a certain sum involving the sine function. For small values of x, we 
know that sin x is approximately equal to x. We need a crude estimate 
between these two values. To get it, recall that for 0 < y S 2/2 we have 


y 
y-3, <siny<y, 
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whence 


a 


2 * 
sin 

ee. <— <1. 

6 y 
The smallest value of the term on the left occurs when y = 2/2, and thus 
we obtain, in our interval, 

y < 2:-sin y. 
Lemma. For any number x (real) not equal to an integer, we have 


|sin tx| R(x) 1— R(x)’ 


and hence for irrational a, 


|sin aoua< nai R(na) 


Proof. The expression |sin zx| depends only on R(x), and hence we 
may assume without loss of generality that 0 < x < 1. Suppose first that 
x = R(x) $4. Let y= 2x < 72/2. Then by the preceding remark, 


SS Ce ee a | 


———_—. = —-— < —- = — < —_., 
lsinzx| |siny| y mx R(x) 


If, on the other hand, 4 < R(x) < 1, then R(—x) = 1 — R(x), and we note 
that |sin xx| = |sin(— oy Thus we can use the same kind of estimate to 
finish the proof of our lemma. 


Theorem 3. Let the hypotheses be as in Theorem 2, and again N 2 Bo. 
Then 


< 4N log N + 40Ng(N) + 2Ko. 
|sin amu 


Proof. This is an immediate consequence of Theorem 2 and the 
lemma. 
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Let F be a polynomial with real coefficients. It is a classical problem to 
estimate sums of type 


N 
a(N, F) = 2, e2tiF(n), 
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Each term in the sum has absolute value 1, and there are N terms in the 
sum. However, depending on the special nature of the coefficients of F, 
one expects a great deal of cancellation to take place in this sum, and 
one expects the sum to have rather small absolute value compared to N. 
We shall first prove a result which suggests fruitful generalizations, and is 
due to Vinogradov (cf. [27], p. 5). 


Theorem 4. Let h be a positive increasing function to infinity. Let k be 
an integer >0. Let 


N 
o(N, x) = Y e2tinkx 
nw 


The set of numbers x with 0 < x <1 such that 
|a(N, x)| 2 N*?H(N), all N > No(x), 
has measure 0. 


Proof. For each positive integer No let Ay, be the set of numbers x 
such that 


|a(N, x)| 2 N*?A(N), all N > No. 


We have 
|a(N, x)|? = o(N, x)o(N, x), 


whence 


. N 1 
| |a(N, x)? ha SS | e2ti(nk—m*)x dx. 
: 0 


m=1 


i= 


But for any integer p # 0 we have 


1 
| a aki), 


0 


whence for N > No, 


1 
N =| |a(N, x)|* dx 
0 
= Nh(N)?: measure(Ay,). 


If measure(Ay,) #0, we obtain a contradiction by taking N sufficiently 
large. Finally, we let A = () Ay, be the union of all Ay, for 


Nor 1 2335s 
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Then A has a measure 0, and A contains the set of numbers x of the 
theorem. 


From Vinogradov’s theorem, we see that given «€ >0, the absolute 
value of the sum a(N, x) is O(N*?**) for almost all x. This leaves open 
the problem of giving a similar estimate for specific numbers x. According 
to our general pattern, such an estimate should arise from a type for x. 
However, only a special case is known, namely when k = 2. To handle 
this case, we shall reduce the sum to the sum discussed in the preceding 
section. 


Lemma. Let «, B be real and « irrational. Let 


F(n) = an? + B, 
and 
N * 
o(N, F) = e2tiF(n) 
n=1 
Then 


N 
N, F)? < 
lo, AP SN+4) naa 


Proof. We have 


st N N 
|o(N, F)|? = 0(N, F)o(N, F) = YY e2tHPoo-Fom), 


n=1 m=1 


In this sum, we consider separately all terms with m =n. Each such term 
is equal to 1, and there are N such terms. The other terms are those 
with n<m and m<n. Thus: 


|a(N, F)|? = N + ys e2 TF) — Fm) 3 e 27i(F(n)— F(m)) 
4 mon An 
(1) | 
SNe RE 3 e2nilF(n)—F(m) 


m<n 


where Re means real part. 
We write n= m-+y. The map 


(n, m) +> (n — m, m) = (v, m) 
gives a bijection of the set of all pairs (n, m) with 


l<m<n<Nn, 


44 ESTIMATES OF AVERAGING SUMS 


and the set of all pairs (v, m) with 
l<v<N and 1<msN-y. 
Furthermore, 
F(n) — F(m) = F(m + v) — F(m) = v?a + 2mva. 


Hence 


(2) ‘iy e2ni(F(n)—F(m)) = ss — e27ilv7a) p2ni2mva 


m<n v=1 m=1 


We shall estimate this last double sum from above by 


where r= e*7"7. 


@ Pe 


(IIT, §3) 


(We note that the term e?"”* has absolute value 1, and thus disappears 


in the estimate.) Now 
Bot Page peer (le ra corto | yen 


and for any real x, we have |sin x| < |1 — e’*|. Hence 


N-v 2: - 9) 


|1 —r| ~ |sin 4xve|’ 


(4) 


rss 


m=1 
Using (1), (2), (3), (4) we see that our lemma is proved. 
Theorem 5, Let «, B be real, « irrational, and let 


F(n) = an? + p. 


Assume that « is of principal cotype <g for all numbers = By. Then 


for N = By we have 


|a(N, F)|? < N + 64N log(4N) + 640Ng(4N) + 4Ko. 


Proof. We apply the lemma, and Theorem 3, with the obvious 


inequality 
4N i] 


|sin ames ~ =i [sin mal’ 
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Remark 1. In practice, the cotype g is a slowly increasing function, so 
that the estimate in Theorem 5 implies in particular 


o(N, F) = O(N ?**) 
whenever g(t) = O(t‘) for some € > 0. 


Remark 2. It would be very desirable to generalize Theorem 5 to 
polynomials of degree > 2. The technique of taking the square is due 
to Weyl [29]. When it is applied repeatedly to handle polynomials of 
higher degree, one obtains a rather poor estimate, far removed from the 
N*?*©« which Theorem 4 leads one to expect. As for the rest of the 
proof, Ostrowski was the first to point out the connection between the 
continued fraction of « and the estimate of the exponential sum [20]. 
His work was continued by Behnke [4], whose ideas we have used, but 
combined with new ones to lead to the very sharp estimates of Theorems 
2 through 5. 


For another approach to the exponential sums, cf. Hardy—Littlewood 
[11], and for other methods, Vinogradov [27]. 

For applications to the next section, we derive one more theorem on 
our exponential sums. 


Theorem 6. Let «, B be real, « irrational. Let m be a positive integer, 
and i 
F,,(n) = man* + B. 


Assume that « is of principal cotype <g for all numbers = By. Then 
for N = Bo, we have 


|a(N, F,,)| < C(N log N)!?(m log m)!? + CN'?m'?g(mN)'? + CKQ?, 
where C is an absolute constant. 


Proof. We apply the lemma and Theorem 3, just as in the proof of 
Theorem 5. 


Remark. The number f is totally irrelevant for the final estimate, 
which is uniformly independent of £. 
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Let h be some function, periodic of period 1, and suppose that h is 
continuous. We can associate with h a Fourier series 


oO 
hoe cae, 
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where c,, is the m-th Fourier coefficient, defined by the scalar product 


1 
a | h(t)e" 27" dt. 
0 


It is proved in any basic course in analysis that if c,,= 0 for all m, then 
h =0 (this being essentially a variant of the Weierstrass approximation 
theorem), and thus the Fourier series determines h uniquely. Let us 
assume in addition that c,,=O(1/m?) for moo (the constant in O 
depending on h). We shall say in that case that the Fourier coefficients 
of h tend to 0 like 1/m?. Then it is an exercise to prove that h(t) is equal 
to the value of its Fourier series for all t (since we assume h continuous). 
If for instance h is twice continuously differentiable, then its Fourier 
coefficients satisfy this condition, as one sees at once by integrating by 
parts twice. 

The estimate of Theorem 6 allows us to treat more general sums, as 
follows. Let F(n) = an? + B be as in Theorem 6. We form the sum 


N N 
S(N, ho F)= ¥ ho F(n)= ¥. h(F()). 
n=1 n=1 
Using the Fourier series for h, we find that 


N © 
SCV, cE) coe” 


00 N 
(1) = » Cm »y e2timF(n) 
00 n=1 


= CaN at 2a Cmo(N, mF), 
m# 


where o(N, mF) is the sum considered in the previous section, and 


oe |, h(t) dt. 


0 


Thus the study of S(N,ho F) is reduced to that of the Fourier coeffi- 
cients c,,, and of our previous sum o(N, mF). 

Observe that if a is of principal cotype <g, then so is —a, because 
the inequality |q« — p| < 1/q can also be written |—ga + p| < t/q. Thus 
in estimating o(N,mF) we may take m to be positive. Under the hy- 
potheses of Theorem 6, and the present assumption that the Fourier 
coefficients of h tend to 0 like 1/m?, we find the estimate 


(oe) 1/2 
|S(N, ho F) ~ coN| S C,(N log N)¥? + C,N¥2 ¥° aun) ES aC Ke 
i J m2 14*0 
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with an absolute constant C,. When g grows slowly, then the error sum 


00 1/2 
zy) = 5 


grows correspondingly slowly, and thus the sum |S(N,ho F)—c)N| 
differs from N*? by a slowly growing function. As an example, we state 
a theorem by selecting special conditions on g. 
Theorem 7. Let h be a continuous function, periodic of period 1, with 
Fourier coefficients tending to 0 like 1/m?. Let co be its 0-th Fourier 


coefficient. Let F(n) = an? + B, with irrational « of principal cotype 
<g for all numbers = By. 


(i) If g(t) = O(t*) for every € >0, t—- then 
|S(N, ho F) —c,N| < KN?+°, 
(ii) If g(t) = O(log t)** for some number k > 0, then 
|S(N, ho F) — coN| S KN"? (log N)!2*#¢, 
(The constant K depends on g, h, and €.) 


Proof. This is an immediate consequence of the estimate for E(N) 
obtained above. The whole point here is that the series 


1 
1 ms 


iMs 


converges for s > 1. 


We observe that from a general point of view, the condition that the 
Fourier coefficients of h tend to 0 like 1/m? is not unnatural: Any 
sufficiently smooth function will have this property, as we have already 
remarked. 

However, in some applications one must deal with functions h which 
are not continuous, and thus whose Fourier coefficients do not converge 
to 0 sufficiently rapidly to apply the previous arguments. One important 
example is the function 


h(t) = R(®), 


where R has been considered before: We have R(t) =t for 0 <t <1, and 
extend R by periodicity to all real numbers. In that case, one tries 
to approximate R by a continuous function whose Fourier coefficients 
behave better. (Cf. Behnke [4], p. 306.) The Fourier coefficients of R 
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tend to 0 like 1/m. If one perturbs R in a 6-neighborhood of 0 to make 
it continuous, and approximates R by a function R,, then the best one 
sees how to do at the moment is to have the Fourier coefficients of R; 
tend to 0 like 1/m?6. By this technique, it is then possible to obtain 
some kind of an estimate for S(N, Ro F) — 3N, namely with an exponent 
3 instead of 4 in situations analogous to those of Theorem 7. We leave it 
as a problem to the reader to find the right method which will yield the 
exponent 4. Since the approximation technique is interesting for its own 
sake, we shall carry out the arguments giving the partial result. 
The graph of R between 0 and 1 looks like this: 


1 


0 1 


For each small positive number 6, we shall define two functions U; and 
L; such that for all t we have 


L(t) = R® S U;(t). 
These functions L; and U; will be continuous, quite close to R except in 


a 6-neighborhood of 0, and will have Fourier coefficients tending to 0 
like 1/m?6. We let U; be the function whose graph looks like this: 


In order words, 


For other 4 U; is defined by periodicity. Then U, is continuous. Fur- 
thermore, if we denote its Fourier coefficients by c,,(5), then by elemen- 
tary integrations, we find 


NI & 


+4 


N| = 


Co(d) = 
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and for m #0, 
1: oe 2nimd a | 


CmlO) = Be ea 


so that in particular, for m 4 0, we have 


i) 
| 
len(9)1 $5 
We have 


S(N, U; 0 F) = co(6)N + ». Cm(O)a(N, mF). 
mst 
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We observe that co(5)N differs from 5N by 6N. On the other hand, by 


Theorem 6 of §3, we obtain an estimate for 
Cn(0)o(N, mF). 
We select 6 = N-"*. We then find that 
|S(N, U;.0 F) — 4N| S N3*E,(N) 


with an error E,(N) similar to E(N). 


One gets a similar estimate with a suitable function L; replacing U;, 


and since 
S(N, L;° F) S$ S(N, Ro F) S S(N, U; © F), 


we see that the estimates for the two extreme sums imply a similar 


estimate for S(N, Ro F). 


CHAPTER IV 


Quadratic Irrationalities 


IV, §1. QUADRATIC NUMBERS AND PERIODICITY 


An algebraic number « is a number which is a root of a polynomial 
equation: 
a," + a,,a"* 4+ +++ ay = 0, 


with rational coefficients a;, not all 0. Multiplying this equation by a 
common denominator for all a;, we see that « is also a root of a 
polynomial equation with integer coefficients, not all 0. In this chapter, 
we are concerned with a special case of algebraic numbers, namely those 
which satisfy a quadratic equation (i.e. those which are a root of a 
polynomial of degree 2). 

Let d be a positive integer. The set of numbers 


x+y/d 


with rational x, y is a field, as is immediately verified. We shall assume 
that it is not equal to Q, and thus that d is not a square. Then without 
loss of generality, we shall always assume that d is square free, i.e. that the 
prime factorization of d does not contain the square of any prime number. 
We note that 1, Jd are linearly independent over Q (otherwise Jd 
would be rational), and thus every element of our field has a unique 
one of the form x + y./d with x, ye€Q. We denote our field by 
Q(/ a). 
Ifa=x+ y/d with x, ye Q, then we define 


a! = x — ys/d, 
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and call « the conjugate of «. The reader will easily verify by direct 
computation that for a, Be Q(</d) we have 


(a+ B) =a’ + p’ and (ob) =2 Pp. 
We define the trace of « to be 
Tr(a) =a+ a’, 


and its norm to be 
N(a) = aa’. 


Both the trace and norm are rational numbers (obvious). We note that a 
is a root of the quadratic polynomial 


(X — a)(X — a’) = X* — Tr(a)X + N(q). 
The quadratic equation for a can also be written in the form 
oo? — 2x + (x? — dy”) = 0. 
We can obtain an equation with integral coefficients, by letting 


x? ye and —2x% =— 
a a 


where a, b, c are integers, relatively prime, and a> 0. Then « is a root of 
the quadratic equation 
aa? + ba+c=0 


with integer coefficients a,b,c, which are relatively prime, a>0, and 
a, b,c are uniquely determined by these conditions. 
We define the discriminant of « to be 


D(a) = b? — 4ac = 4a’y7d. 
Since « is real, irrational, we have D(a) > 0. 
We shall say that « is reduced if «>1 and —1 <a’ <0 (or equiva- 


lently, —1/a’ > 1). If « is reduced, so is —1/a’. 


Theorem 1. Given a positive integer D, there is only a finite number of 
reduced elements of Q(./d) whose discriminant is D. 
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Proof. Let « be reduced, and have D as discriminant. We have 


= —b— D 
oe and “| ee no. 
2a 2a 


where € = 1 or € = —1. If €e = —1, then 
—-b—-./D>0 and —b+./D<0, 


a contradiction. Hence e = 1. Then from these same conditions, 


—b+./D>2a>b+./D. 


Hence b <0 and 0< —b< ED: Hence |b| is bounded in terms of D, 
and the preceding inequality shows that a is bounded in terms of D also. 
Thus finally |c| is bounded, and our theorem is proved. 


Lemma 1. If « has the discriminant D> 0, and B is equivalent to a, 
then B has the same discriminant D. If m is an integer, and B =a +m, 
then again B has the same discriminant as «. 


Proof. Direct brute force computations. 
Theorem 2. Let « be real, quadratic, irrational. 
(i) In the continued fraction 
a= [ao, Bleue 2An—1 ? oa 


the number a, has the same discriminant as a for n= 1. 

(ii) If « is reduced, then a, is reduced for all n= 1. 

(ii) If a is not necessarily reduced, then a, is reduced for all n 
sufficiently large. 


Proof. From the lemma, we conclude at once that «, has the same 
discriminant as « (using the construction of the continued fraction, and 
Chapter I, §3). This proves (i). Furthermore, «, > 1. If a is reduced, and 


1 
a=at— 


with an integer a, and 6 > 1, then —1/f’ =a— a’ >1, so that B is also 
reduced. This proves (ii). To prove (iii), note that 


a’ i Pn-1 
a! = _In-1 wee ire Qn-1 
Gn k is Pr Qn a = Pn 
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For n large, the fractions p,_,/q,-; and p,/q, are close to a, and hence 
a’ — p,-1/q,-, and «' — p,/q, are close to « — a Hence the expression on 
the right is negative. On the other hand, 


is small, and hence we see that a, + 1 >0. This proves (iii). 


Let « be any (real) irrational number. We shall say that its continued 
fraction 


lagnayaee. | 


is periodic if there exists an integer k such that for all n sufficiently large, 
n+, = 4,. We shall say that « is purely periodic if a,,, =a, for all n, and 
we say that k is the primitive period if it is the smallest positive integer 
with this property. 

We shall use the following standard notation for periodic continued 
fractions. We write 


[ane 4, Geen s+ Gaga 
as an abbreviation for 
[ao, vp, Appts oes sOptis Aptis +++ srt, neal 


with a repetition of a,,,,...,4,4,, so that for n 2r+ 1 we have a,,, =a,. 
A continued fraction is purely periodic if and only if, in the preceding 
notation, it can be written in the form 


cred 


Lemma 2. Let « be reduced, and let a be an integer. Let 


a=at+—. 
ay 


Then «, is reduced if and only if a<a<a+1, that is, a is the 
greatest integer S a. 
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Proof. If a satisfies the stated condition a<a<a+t 1; then ay> 1, 
and a’, = 1/(a’ — a), so a, is reduced since a2 1 and « <0. Conversely, 
if «<a then «, <0 and ifa+1<a then a, <1, so that a, cannot be 
reduced. 


We observe that the relationship between a and a, in Lemma 2 
determines each number from the other. Indeed, we have 


Thus, we see that a is also the greatest integer < —1/ax,. Furthermore, 
aw’, and hence «, is uniquely determined by «,, and hence also by «,. 


Theorem 3 (Euler-Lagrange). Let « be real irrational. The continued 
fraction of « is periodic if and only if « is quadratic. Assume this is the 
case. Then « is reduced if and only if its continued fraction is purely 
periodic. 


Proof. Assume that « is quadratic irrational. We know from Theorem 
2 that a, is reduced for large n, and by Theorem 1, combined with 
Lemma 1, there is only a finite number of possible values for such «,. 
Hence for some n and k > 1, we have a, = «,4,. From this, the periodic- 
ity follows at once. Suppose in addition that « is reduced, and that 
Ln = m+, for some m>0, kK21. From Lemma 2, and the remarks 
following Lemma 2, we conclude that «,,_, 1s uniquely determined by «,,, 
whence also &,-; =Om+,-1- From the definition of the continued frac- 
tion, we conclude that it is purely periodic. 

Conversely, if the continued fraction is purely periodic, we can write 


Cage eae | dose. «Aaa Gol: 
From the relation «=o, we see at once that « is quadratic, and 


furthermore, « =a, for arbitrarily large n. By Theorem 2(iii) it follows 
that « is reduced. If the continued fraction is merely periodic, with say 


Ci CARRERE ey rcneiartita |: 


then 


a+ = (4,41, --29teee 


is purely periodic, whence quadratic. Since «,,, is equivalent to «a, it 
follows that « is quadratic. 


fLV582) UNITS AND CONTINUED FRACTIONS 55 


EXERCISE 


Show that if a is reduced, and 


“g£= [@o; eee ays I 
then 


i a 
ee, = las; db 506 1- 
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Let a be quadratic as before. We shall say that a is an algebraic integer 
if Tr(a) and N(«) are integers. If a is an algebraic integer, then its 
quadratic equation 


X? — Tr(a)X + N(a) =0 


has integer coefficients, and conversely. We shall write « in the form 


utv/d 
a 


with u, v€ Q. The reason for this is that we have Tr(a) = u, and 


u* — dv* 
N@ ==>. 


Theorem 4. An element « of Q(./d) is an algebraic integer if and only 
if u and v are integers, and the following conditions are satisfied: 


‘ = v (mod 2) if d = 1 (mod 4), 
(+) u=v=0(mod2) if d=2,3(mod 4). 

Proof. If the stated conditions are satisfied, then one verifies at once 
from the definitions that « is an algebraic integer. Conversely, if a is an 
algebraic integer, then from its trace we see that u must be an integer. 
Then 

—dv* = 4N(a) — uv? 


is also an integer, and since d is assumed square-free, it follows that v is 
also an integer. Then 


u? — dv? =0 (mod 4). 


56 QUADRATIC IRRATIONALITIES [IV, §2] 


From this congruence, we deduce at once that the congruence conditions 
on u,v stated in the theorem must hold. 


Let 
d= 


tld if d= 1(mod 4), 


0=/d if d = 2,3 (mod 4). 
Then from Theorem 4, we obtain: 


Corollary 1. An element « of Q(/d) is an algebraic integer if and only 
if « can be written in the form 


a=x+y0 
with integers x, y. 
Corollary 2. The algebraic integers in Q(./d) form a ring. 


Proof. From the representation in Corollary 1, it suffices to prove that 
6? is an algebraic integer, and this is clear. 


A unit in Q(./d) is an algebraic integer whose inverse is also an 
algebraic integer. If w is a unit, then from the relation ww™' = 1, we 
conclude by taking the norm that N(w) = +1. Conversely, if aw’ = +1, 
then @ is a unit, and @' = +q@’. Thus the units form the multiplicative 
group of non-zero algebraic integers whose norm is 1 or —1, ie. form 
the set of solutions of the equation (called Pell’s equation) 


u? — dv? 

—_—— = +l, 
4 ats 

where u, v are integers, satisfying conditions (*) of Theorem 4. 
Let U be the group of units. We map U into R? by 


L: @ + (log|a|, log|q’|). 


Then L is a homomorphism of U into the straight line satisfying the 
equation X + Y=0. The kernel of L consists of those units having 
absolute value 1, and the quadratic equations for such units have integer 
coefficients of absolute value < 2. Hence there is only a finite number of 
elements in the kernel of L, which is therefore a finite group, i.e. a group 
of roots of unity. Since the only roots of unity which are real are +1, it 
follows that the kernel of L is {1,1}. Furthermore, in any bounded 
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region of R? there is only a finite number of elements of the image of L, 
because any bound on log|@| and log|w’| yields a corresponding bound 
for the integer coefficients of the quadratic equation satisfied by w. We 
shall see in a moment that the image of L contains vectors other than 
the zero vector. 


Let w, be a unit #1, —1 such that L(w,) is a vector of smallest length 
in the image of L. Then.we contend that every element of the image of 
L is an integral multiple of L(w,). 


Proof. Let W, = L(q@,), and let W be in the image of L. Then there is 
a real number ft such that W=tW,. Let m be an integer such that 
ms<xt<m+41. Then 


W—mW, =(t—m)W,, 


and the length of W— mW, is smaller than the length of W,. Hence 
WwW — mW, = 0, so W =mW,, as was to be shown. 


When we have shown that there exist non-trivial units (i.e. units other 
than 1 or —1), we therefore obtain the following theorem: 


Theorem 5. Modulo {1, —1}, the group of units is infinite cyclic. There 
exists a unique unit w, which generates this group and such that w, > 1. 


Proof. A unit w, such that L(w,) is a generator for the image of L is 
of course not uniquely determined, and any one of the units 


al al 
M1; —@M,,@, 3 aC 


would also have this property. However, exactly one among these four 
numbers will be > 1, thus proving our theorem. 


We shall now determine non-trivial units. We let D be the discriminant 
of the number @ defined previously, so that by a direct computation, we 
find 

a d if d=1 (mod 4), 
~ (4d if d= 2,3 (mod 4). 


We shall call D the discriminant of the field Q(/d). From Theorem 4, 
we obtain: 
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Theorem 6. An element a of Q(./d) is an algebraic integer if and only 
if it can be written in the form 


u+v./D 
> an 


with integers u,v such that u= Dv (mod 2). Furthermore, the algebraic 
integer a is a unit if and only if, in this representation, we have 


u* — Dv? 


= +1. 
; + 


From the characterization of units in Theorem 6, we shall deduce 
results of Pell, Euler, and Lagrange. 


Theorem 7. Let « be a reduced number in (./d), having D as dis- 
criminant [where D is the discriminant of Q(./d)]}. Let k be a period of 
its continued fraction, let v be the greatest common divisor of 


(Qi—1> Pe-1 — V-2> Pr-2) 
and let 


U = Dy-1 + Gx-2- 


Fp ee 


Then 


is a unit > 1, and every unit > 1 is of that type. We have 
N(@) = (— 1. 

Proof. In the notation of Chapter I, §3, we have 
= O,10 


since by definition 


Py-1% + Dy-2 


® =[Go,....,0.4e— . 
Ge-1% + Ye-2 


Thus « is a root of the equation 


Qk-1 a? + (Qe-2 — Pr-1)% — De-2 = O, 


and, as before, we let 


aw? + ba+c=0 
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be the equation for a with relatively prime coefficients a,b,c and a> 0. 
Then we can write 


(1) Qk-1 = av, Qx-2 — Py-1 = bv, i — Px-2 = CV. 


As in the theorem, let 


(2) Uppy + 94-2- 
Then 
u — bv 
Py-1 = ) and Q-1 = av, 
(3) 
u+ bv 
Py-2 = — Cv and Q-2 = ay 


From the relation 
G1 Pe—-2 — Pr-1 9-2 = (— 1D eae 


we conclude: that 
u* — Dv? = (— 1)*4, 


whence the norm of @ is equal to (—1)*. Since the trace of w is equal to 
u, we conclude that @ is a unit, w > 1. 


Conversely, let 
u+v,/D 
(60) => =i 


be a positive unit > 1, with integers u,v. Then u,v 2 1 (because among 
the four numbers 


exactly one is > 1). Let p/q and p’/q’ be defined in a way similar to (3), 
namely, 

u— bv 

ae 


q = av, 


; ,  u+ bv 
p’ = —cv, rca a 


Since’ u=vD (mod 2) by Theorem 6, and b=D (mod 2) because 
D? = b? — 4ac, it follows that p, p’, q, q’ are integers. Then 


Gq pap 
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are proportional to a, b,c and hence « satisfies the equation 
qu” + (q’ — p)a + (—p') = 9, 


or in other words, 
pa + p’ 
oh -. 
qa + q 


We have trivially, using the definition D = b? — 4ac, 
pq’ — gp’ = N(w) = +1. 
Using the inequalities of Theorem 1, §1, we find: 
, ut bv u —v./D , N@) 
2 Z oD) 
0 if NW) =1, 
—-1 if N@=-—1. 


, @a—bjv—u v/D—u_ ,_ ~N(@) 
=! ie 5} = 5} = = 
—-1 if N(@@)=1, 
0 if N(@) = —-1. 


q 


Case 1. N(@) = —1. Then g>q' 20. If q’ =0, then p’ = 1 and 
“= i 
=p a 


so that a =a, and we are in the situation of the first part of the proof, 
with k= 1. If q’ > 0, then we can apply Theorem 7 of Chapter I, §3, to 
conclude again that « = o,_,a, and that we are in the same situation as 
the first part of the proof. 


Case 2. N(w)= 1. Then g’>0 and q2q’. If q=q' then q=q' =1 
because q, q’ are relatively prime. Then a trivial computation yields 


1 
1+ 1/a 


v 


a= p+ 


so that we are again in the situation of the first part of the proof, namely 


a=[p’,1,a]. If q>q’, then we can apply Theorem 7 of Chapter I, §3, 
once more to achieve the same result. 
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We have shown that in all cases, w is constructed from « as in the 
first part of the proof, which therefore yields all units > 1. 


We shall obtain finally a slightly different description of the units 
constructed in Theorem 7. 


Theorem 8. Let a be reduced, with discriminant D [the discriminant of 
Q(/d)]. Let k be a primitive period. For each integer m= there is a 


unique unit @,, Such that - 
e 4 1 
= =O 


(where o,, is the n-th continued transformation of a). We have o,, = 7, 
and w, > 1 generates the group of positive units. 


Proof. From the definitions, we see that w,, is precisely the unit 
constructed in the first part of the proof of Theorem 7, and by Chapter I, 
§3, we have 


pean m 
Omk-1 = 9x-1> 
whence w,, = w7. This proves our theorem. 


We see that the positive units > 1 can be interpreted as eigenvalues 
of certain linear maps, namely the continued transformations of « cor- 
responding to the periods. 

To apply Theorem 7 or 8, one must start with a reduced number «. 
The next result gives a simple way of finding such a number. 


Theorem 9. Let 0 be the number defined for Corollary 1 of Theorem 4. 
Let [6] be the greatest integer < 6. Then 


1 


=u 


is reduced. 


Proof. The proof will be left as an exercise. 


IV, §3. THE BASIC ASYMPTOTIC ESTIMATE 
We have seen in Chapter II, §2, that quadratic numbers are of constant 


type. Thus the asymptotic estimate of Chapter II, §3, applies, except for 
the special case of the inequality 


c 
0<qa—p<- 
q 
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with a constant c. We shall take care of this case by using the special 
properties of quadratic numbers. We follow [15]. 
Theorem 10. Let « be a real quadratic irrational number. Let c be a 
positive number such that the inequality 


c 
O0<qa—p<- 
q 


has infinitely many solutions in integers q, p, and q>0. Let A(N) be the 
number of solutions of this inequality for 1<q<N. Then there exist 
numbers c,,C, > such that for all N we have 


|A(N) — c, log N| S cp. 
In other words, 
A(N) = c, log N + O(1). 


The proof will need some lemmas. 
Let ae Q(/d) with d > 0, square-free as usual. We let as before 


d if d = 2,3 (mod 4), 


@= 
: ul if d= 1 (mod 4). 


Then 1, 6 form a basis of the algebraic integers of Q(/d) over Z. We 
let 


a=a0+b 


with rational a,b and a #0. Let e be a positive integer such that ea and 
eb are both integers. Then ex is an algebraic integer. 

We shall use the phrase “sufficiently large (resp. small)” to mean 
“greater (resp. smaller) than a constant depending only on 06, « and c”. 


Lemma 1. There exists an integer k > 0 having the following property. 
An integer q (sufficiently large) is such that 


(*) O<qu— p<? 


for some integer p if and only if there exists an integer p such that 
qu — p is positive, sufficiently small, and 


k 
e 


(+x) IN(q« — p)| S 
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Proof. We distinguish cases, depending on whether ce?(«’ — a) is or is 
not an integer. 
Suppose that ce?(a’ — a) is not an integer. Then we take 
k = [ce*(a’ — a)]. 


First assume (x). Then 


: 2 
ce 
IN(qea — ep)| < eas — p| = ce? 


a” — 4 
q 
If p/q is close to a, then a — p/q is close to «’ — a. Since the norm of an 
algebraic integer is an integer, we conclude that 
|N(qea — ep)| Sk, 
thereby proving (xx). 


Secondly, assume (xx), and also that qa —p is positive, sufficiently 
small. Then 


The quotient k/ce?|«’ — «| is a fixed number < 1. For q sufficiently large, 
a’ — p/q is close to a — a, and hence the right-hand side of our inequality 
is < c/q, thereby proving (x). 

Suppose that ce?(«’ — a) is an integer. We distinguish two subcases, 
depending on whether a < a’ or a’ <a. 

If « < a’, then for q sufficiently large, we have 


Lee 
q 


We take k as before, and the first part of the argument runs as before, 
because 


IN(gea — ep)| < ce? (« 2 r). 
For the second part, we now use the fact that 


ee ey 
q 
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so that 
0 ee ee 
Se PS Tce —a) q’ 


If « <a, then for q sufficiently large, we have 


Cea, 


This time, we take 
k = ce?|a’ —a| —1. 


In the first part of the argument we have 


Pp 
v | <|a’—al, 


and the desired conclusion follows. The second part of the argument is 
carried out as before, thereby proving the lemma. 


Remark. Any c 21 in Theorem 10 will do. Furthermore, Lemma 1 
and the fact that the norm of an algebraic integer must be an integer 
show precisely how small we can take c and still get infinitely many 
solutions. 


In view of Lemma 1, we must count the number of integers q sat- 
isfying 1<q<N, such that there exists p for which qa — p is positive, 
sufficiently small, and 


|N(qea — ep)| Sk. 


Let m be an integer, 1 <m<k. We shall prove that our asymptotic 
estimate holds for the number of solutions of 


|N(gea — ep)| =m, lsqgsN, 
with ga — p positive, sufficiently small, provided that there is at least one 
solution. Adding up these estimates for m=1,...,k and using the fact 
that our original inequality actually has infinitely many solutions, we 
obviously obtain a proof of our theorem. 

We shall see that our problem can be reduced to counting certain 
units. For the rest of the proof, it is convenient to define a certain 
equivalence relation. Let &,7 be algebraic integers. We shall say that 
they are equivalent if there exists a positive unit @ such that ce =n 
This is obviously an equivalence relation. 
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Lemma 2. Given a number B > 0, there exists only a finite number of 
inequivalent algebraic integers € such that |N(é)| < B. 


Proof. Let wp be a generator > 1 for the positive units. If |éé’| < B, 
we shall see that there is some power w), such that 


|ao¢| and —|(@G 


are bounded in terms of-B. Since there is only a finite number of 
algebraic integers with bounded conjugates, this will prove what we want. 
Suppose that one of the two absolute values |€|, |¢’| is > @ 9. By symme- 
try, say |€| > @o. Let n be the integer such that 


Gn icl= Os. 
Then 
[¢/a@5| < ao. 


On the other hand, |¢’| < B/|é|, and 1/|€| < 1/w§, whence 


iS Bo. 
of] ~ laos 
This proves what we wanted (with r = —n). 


Lemma 3. Let qo, Po be integers, qo #0, and let €) = qoex — ep. 
The set of units u such that ué) can be written in the form qea — ep 
with integers q, p is a group. 


Proof. Write a, = ea and b, =eb. We write a unit u as u=x0+y, 
with integers x, y. We shall prove that the condition stated in the lemma 
is equivalent with a congruence condition on x. 

We have 

Co = 40419 + Gob; — po. 
Then 


UE = (x(qoby — ePo) + ¥4041)% + Xqoaid + Y(dob: — ePo), 
and we must find a necessary and sufficient condition that this expression 
is of type 
qa,a + gb, — ep, 
with integers q, p. This amounts to the pair of conditions 
X(dob, — €Po) + Y4041 = 941; 
XGo4,d + y(qob, — ePo) = 9b; — ep. 


66 QUADRATIC IRRATIONALITIES [1V"35) 


The first one simply means that a, divides the left-hand side. Let w be 
the g.c.d. of a, and (qob, — epo). Let a; = way. Then the first condition 
is equivalent with ay|x (provided qob, — epo #0, a case we leave to the 
reader). We shall write x = agx*. 

The first condition being satisfied, our second condition yields another 
divisibility condition, namely that e divides 


— 
x (<odoaid —b, aah — 20), 


Let t be the g.c.d of e and the expression in parentheses which we have 
just obtained. Write e = teo. Then our last condition amounts to 


@g|X*. 
Hence finally, our two conditions are equivalent with the divisibility 
Ayeo|x. 


Since u=! = +u’, it now follows at once that the units satisfying our 
divisibility condition form a group, as contended. 


Lemma 4. In Lemma 3, assume that €, >0, and let S be the subgroup 
of positive units in Lemma 3. Assume that S is infinite cyclic, and let w 
be a generator of S,0<a@<1. Write a unit ue S in the form u=@". 
The subset of S consisting of those units u such that u is sufficiently 
small, and such that in the expression u€y = qea — ep we have q > 0, is 
of one of the following types: It is empty, it consists of all sufficiently 
large n, or all sufficiently large even n, or all sufficiently large odd n. 


Proof. We have 
Udo — U'lo = ge(a — a’). 


When uw is small, then ug, is also small, and u’ is large (positive or 
negative) because N(u) = uu’ = +1. Hence —u’€5 is of the same order of 
magnitude as ge(a — a’). Suppose first that a>’. Then g>O0 if and 
only if —u’E >0 whenever u is sufficiently small. Let « = wa’, and let 
n = —o. The condition —u’&) > 0 is then equivalent with the condition 
€"7 > 0 (write u’ = q@™ and multiply the inequality by the positive num- 
ber "). We then have four cases to consider, depending on whether 
e=1, €=—1, 7 >0, or 7<0. These four cases obviously give rise to 
the four possibilities described in the lemma. We note finally that the 


argument is entirely similar in the case that «’ >, thus proving our 
lemma. 
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Lemma 5. Let qo > 0 and py be integers, and let 
So = oed — ePo 
be >0. Let Ao(N) be the number of pairs of integers p, q with 
beq=N 


such that, if we let € = qea — ep, then: 


(1) € is equivalent to 9; 
(2) &€>0; and 
(3) € is sufficiently small. 


Then there exists a constant co = 0 such that 1)(N) = co log N + O(1). 


Proof. Let S be the subgroup of positive units in Lemma 3. If S has 
only one element, we can take cy =0, and A,(N) is then bounded, i.e. 
O(1). Assume that S has more than one element. Then S is infinite 
cyclic. Let w be a generator with O<@<1. An integer g such that 
€ = qea — ep satisfies the three conditions of the lemma is then deter- 
mined by €, which can be written in the form € = ué9, with u = w”, and 
we deal with one of the four possible types described in Lemma 4. If the 
set is empty, we are again done. Suppose this is not the case. Consider 
for definiteness the case when we deal with all sufficiently large even 
integers n. From the equation 


ulo — u'So = ge(a — a’), 
we see that there exist two constants k,,k, > 0 such that 
k,qslu'| sk, 
or in other words, 


kiqsu" Skq. 


Given a number k, > 0, we note that the number of units u = @” such 
that n is a positive even integer and 


1 <u =(1/o)! <k,N 


is equal to the number of positive even integers n such that 


log(k3N) _ 1 
(ere) «eat + oO 


and is therefore given by 


Co log N + O(1), 
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where cy = } log(1/m). [We see that the constant k, is absorbed in the 
error term O(1).] 


For positive numbers B, let us denote by S(B) the set of units 
u=q"eS such that n is sufficiently large, and u-' < B. Then for suit- 
able constants k,,k, > 0 we have 


card S(k,N) S A)(N) S card S(k,N). 
Using our preceding estimate, we see that the theorem is proved. 


We note that even though an asymptotic estimate was known for 
almost all numbers, the estimate of Theorem 10 was the first to give a 
similar result for specific numbers which could be exhibited explicitly. 

It is a fundamental problem in the theory of algebraic numbers to 
extend the results of this chapter to those numbers of degree > 2. As 
Liouville noticed, if « is an algebraic number of degree n, then one 
obtains a trivial generalization of the fact that quadratic numbers are of 
constant type, namely there exists a constant c > 0 such that 


G 
|qa — p| a 


The proof is essentially the same as that given for the example of Chap- 
ter II, §2. The best-known type for algebraic numbers is given by the 
Thue—Siegel-Roth theorem [23], but in spite of the difficulties which 
have been encountered historically to reach this result, one expects that 
the type O(t‘) can be replaced by O(log t)’ for some p> 0, and thus 
Roth’s theorem still appears as rather far removed from a good descrip- 
tion of the situation. 

For other possible directions in generalizing diophantine approxima- 
tions to algebraic numbers, cf. Schmidt [26] (for a Liouville type result), 
and also Bernstein [5], [6], [7], and Hasse—Bernstein [12]. 


CHAPTER V 


The Exponential Function 


V, §1. SOME CONTINUED FUNCTIONS 


It is a general problem to determine the continued fractions for values of 
classical functions suitably normalized. We shall describe a solution of 
this problem in a very special case which will allow us in particular to 
get the continued fraction for e. 

We start with the function 


1 aan ieee xe 
GERI ae + 2) 31 


1 
fe,x)=14+—x + +: 


n 


= 1 bd 
7 c(c + 1)--\(c +n—1)n! 


The number c is taken for the moment to be any real number not equal 
to any integer 0, —1, —2,... so that the series is defined. Then it is 
trivial to verify that 


f(c, x) = f(c + 1, x) + ——— fc + 2, x), 


ne 1) 
so that 


fle + 1, x) _ oe 


i ye e+ 20 


er) 
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and hence 
fle+1,x) _ ! 


f(c, x) ry x f(c+2,x) 
* ce + fe t+ 1x) 


This almost looks like a continued fraction, but the 1 is on the wrong 
side in the denominator. Thus we transform this expression by letting 
x = z7, and write 

z f(c + 1,27) 1 


e f(z) ~¢ 2 Ale + Dez.) 


z e+1f(c+1,z7) 


This is now in the form where we can express it in terms of the formal- 
ism of Chapter I, §1, in a convenient way, namely, by induction: 


z f(c+1,2°) _ ie Cae 
c. flc.z2) _ f(c, 27) boa) Z a 9 “at2 | 


where 


? Detn+ lyfe re iz) 
ta” g”StéOf (cs SEED 


We shall therefore obtain the continued fraction expansion for the value 
of the left-hand side if we can substitute special values for c and z such 
that (c + n)/z is an integer 2 1, and such that the last term «,,, is 2 1, 
for all integers n 20. For instance, it is immediately verified that the 
values 


for any integer y = 1 will fulfill these conditions. In this way, we obtain 
the Lambert continued fraction (which actually was already known to 
Euler). For these values of c and z, the function can be expressed in a 
more familiar way. Indeed, we have 
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because the coefficient of (w?)" in both power series is 


1 ee! 
23 + 1)---(G+n—1)4n! — 3-5---(2n + 1)2-4--- In 


a 
ers hl 


2 
cveror(h) 


We obtain the Lambert continued fraction if we let w = 1/y: 


Similarly, 


Theorem 1. For every integer y = 1, we have 


elly ae ety 
a 4 ed LOys; 3y; SVs | 


and, in particular, for y = 2, 


e—l 
an =f. 0,10, 22.1. 


The continued fraction for e will be obtained easily in the next section 
from the continued fraction for (e — 1)/(e + 1). Observe that these two 
numbers are not equivalent, however, since the determinant of 


1 -1 
1 1 
is equal to 2. 


We note that the same recursive method which we used for the expo- 
nential function also works for a wider class of functions (hypergeometric 
functions, as they are called), including the Bessel function. However, 
just as in Theorem 1, one obtains continued fraction representations only 
for special values of the variable, and it is a problem to find the right 
approach to get the answer in general. There is also the difficulty that 
for functions satisfying, say, a second-order linear differential equation, 
like the Bessel function J, one gets information on J’/J, but one still does 
not have much information concerning the continued fraction of J itself 
(for rational or integral variables). 
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V, §2. THE CONTINUED FRACTION FOR e 


We shall give Euler’s argument to prove: 
Theorem 2. The continued fraction of e is 
c= 24a dbl t,.....]. 
In other words, ag = 2, and for m 2 1, 
ds = Gsaene and Asn—-1 = 2M. 
Proof. Let r,/s, be the principal convergents to the number 


eral 


a =| 


Inverting the continued fraction of Theorem 1, we know that 
ial 2. 6910) oe]; 


and obviously, 


eatt! 
Sigal 


Let € = [2, 1, 2, 1, 1,4, 1, 1, 6, Lode...) pemdulet pela, beathesprincipal’con- 
vergents to €. We shall prove that for all n => 0 we have 


(x) Z P3ntt =" t+ Sp and Q3nti = Th — Sp- 


These relations are verified by a direct computation for n=0,1. By 
definition, we have for n = 2: 
r, = (2 + 4n)r,-1; + 1,-2; 
CS al (2 a 4n)Sy—1 + Sn-2- 
On the other hand, we multiply the recursion formulas for p,,q, by the 
number indicated on the right in the following table, and add. 
P3n-3 = P3n-4 + P3n-s 1, 
P3n—2 = P3n—3 + P3n-4 = 
P3n-1 = 2ND3,-2 + Pan-3 2, 
P3n = P3n-1 + P3n-2 1, 
P3nt1 = P3n + Pan-1 1. 
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We obtain (for n 2 2) formulas similar to those for r,, s,, namely 
Panti = (2 + 4n)p3,-2 + D3n-s> 
Gantt = (2 + 4n)3n-2 + Fan-s- 


Relations (x) now follow by induction. We conclude that 


We know that p3,4:/d3n+1 approaches € as n tends to infinity. Since r,/s, 
approaches a, it follows that € = e, as was to be shown. 


By a similar method, one can get the continued fraction for e?” with 
an integer y 21. Cf. Perron [21]. 
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The asymptotic result of Chapter II, §3, is valid only under certain 
conditions determined by the type of the number. We shall now repro- 
duce a result of Adams [1], determining the basic asymptotic estimate for 
e. We let A be the function Az iq considered in Chapter II, and we let 
Ao(N) be the number of solutions in relatively prime integers p,q of the 


inequalities 
(1) 0<qe—p<I1/q and besa N. 
To begin with, we can sharpen Theorem 10 of Chapter I, §4. 


Theorem 3. Any solution p/q in relatively prime integers with q > 0 of 
the inequality |qe — p| < 1/q is a principal convergent. 


Proof. We must show that a solution of the form 


Pn + TPnt+1 
Gat 1 Gass 


with r=1 or r=a,,,—1 cannot occur. If n= 3m—2 or n=3m—1, 
then a,,, = 1 and so r=0. Suppose n = 3m. By the lemma of Theorem 
10, Chapter I, §4, we must show that the only solution of 


C3m+2 — 7 1 


ea 
Gam + €3m+293m+1 3am + 19 3m41 
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is r=0. That r #1 follows from 
Csm+2 > Asm+2 = 2(M+ 1) and 43m+1/43m = UL, 1, 2m,...] < 2. 


One sees that r# 43,4, —1 in a similar way, thereby proving Theo- 
rem 3. 


In what follows, we shall meet the function 4*I'(x + 3), which is 
strictly increasing for x >0. We shall denote its inverse function by g. 
Simple estimates show that g(x) is asymptotic to 


(log x)/log log x. 
Theorem 4. We have 


Ao(N) = 3g(N) + O(1) 
and 
A(N) = §(2g(N))*? + O(g(N)). 


Proof. We shall first obtain a special case of Theorem 4, for special 
values of N. 


Lemma 1. We have for all n, m2 1, 


Ao(Gn) = 20 + O(1), 
A(G3m+1) = (2m)*” + O(m). 


Proof. The first assertion is a direct consequence of Theorem 3, the 
definitions, and the fact that the principal convergents alternate on both 
sides of e by Theorem 5 of Chapter I, §2. (This gives rise to the factor 
3 in front of n.) As for the second assertion, we must determine the 
multiples of the convergents p,,q, which are also solutions of (1), and 
0<q,« — p,. Using once more the Lemma of Theorem 10, Chapter I, 
§4, with r = 0, we see that q = kq, and p = kp, is a solution of (1) if and 
only if 

oe 7 oe 
Qn ate Cn+29n+1 k Qn Cn+2 Qn 


If n= 3m — 1 or n= 3m, then the condition implies k* < 4, so k = 1 is 
the only possibility. If n = 3m — 2, the condition amounts to 


k?<2m+O(1), ie. 15k <(2m)!? + O(1). 


For such k, we note that kq3,,-2 < q3m+1 (for m sufficiently large). Hence 


[V, §3] THE BASIC ASYMPTOTIC ESTIMATE 1S 


modulo O(m), we find 


m-1 


m/2 
Adsnes)= Y Q2)'? = | aga de=sem. 
v= (0) 


(We restricted the sum to v even because we consider only those n such 
that 0 < q,e — p,, and those are the even n, so that here this amounts to 
even v.) Our lemma is proved. 


Proving Theorem 3 now essentially amounts to obtaining m as a 
function of q3,,-2- 


Lemma 2. There exist constants c,,c, > 0 such that 
c,4"T(m + 3) S dame1 S C24" (m + 9). 
Proof. We note that the equations 


Qam+2 = 2(m + 1)dam+i + I3m> 
— W3m+1 = 13m + I3m-1> 


Gam = W3m-1 + I3m-2> 


can be solved to yield 


A3m+1 = 2(2m ie 1) es am-5 


I3m-2 I3m-2 
so that 
dam+t = [2(2m + 1), 2(2m — 1), ...]- 
G3m-2 
Thus 


Iame1 = 2(2m + 1)43m-2 2 27(2m + 1)(2m — 1d3m-s 2°" 


and hence 
Geet l (m+ 3) 


is clear. Conversely, 


| 
W3m+1 xe 2 2 a6 1 oe 
pe = a 1) 


1 
< 2(2m + (1 * (am + 1m — 5) 
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and proceeding inductively we see that 


m 


1 
damsi S 2"(2m + 1)(2m — 1): |] (: * Ay + HQ — 7, 


v=] 


sO 
Fam41 SC24"T(m + >), 
where c, is determined by the infinite product. 


To prove the theorem, we find to any given N the integer m such that 
W3m-2 = No< Oat Thus 


c,4" "Tm —14+3) SN <c,4"T(m + 3), 


and consequently 
g(N/c2)< mS g(N/c,) + 1. 


Since g(x) grows like (log x)/log log x, we conclude that 
m = g(N) + O(1), 
whence Theorem 4 follows from Lemma 1. 


Remark. The procedure here differs only very slightly from that in 
Chapter II, §4, and we could simply have applied what we did before. 
We chose to repeat the arguments ab ovo, as in [1], partly to show how 
one obtains a slightly more accurate bound for gq, (n= 3m-+ 1) in 
Lemma 2 than by the procedure of Chapter II, which would have intro- 
duced an extra factor 2”. We also remark that we could use the function 
x* instead of I(x), since factors c* (with c constant > 0) do not affect the 
asymptotic behavior of the inverse function of c*x*. 


From the continued fraction of e we can also determine a type. 


Theorem 5. If g(x) is as above, the inverse function of 4*T(x + 3), then 
e is of type < 2g + O(1). 


Proof. Given N we find q, such that 


Qn-1 < N = Qn: 
We have 
Aig 4 ee 1. 
Qn-1 n-1 
Hence 
N In 


< <q,_ 
a, #1 pay 
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If n= 3m or 3m — 2, then a,=1. If n= 3m —1, then from the expres- 
sion m = g(N) + O(1) obtained previously, we find 


An = Azm—1 = 2m = 2g(N) + O(1). 


This proves what we wanted. 

We note that the type of Theorem 5 is essentially best possible. 

In [2], Adams has extended the estimate of Theorem 4. He treats the 
case of more general functions W(t) = w(t)/t as described in Chapter II, in 
the range where the result of Theorem 8, Chapter II, §3, does not apply. 
He also treats a wider class of numbers, those whose continued fractions 
look like that of e (including Hurwitz continued fractions—cf. Perron 
[21]). 

The general problem is to investigate values of the classical functions, 
suitably normalized, from the present point of view. (Cf. [17] for a 
general discussion.) The function e‘ is of course the simplest one. The 
first problem which arises, and probably the simplest, is to determine 
the continued fraction for e* where a is rational (or even an arbitrary 
integer). One wants to see how the special analytic properties of one of 
the classical functions are reflected in the arithmetical properties of their 
values. Theorem 1 of §1 gives an example of this reflection, showing that 
all numbers obtained as values of a function suitably normalized, and 
over a suitable domain of definition, have similar continued fractions. 
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APPENDIX A 


Some Computations in 
Diophantine Approximations! 


By W. ADaAMs at Berkeley and S. LANG at New York 


Let w,,...,W,, be (real) numbers, linearly independent over the rationals. 
Let B be an integer >0 and c a number >0. We define A(B,c) to be 
the number of solutions of the inequality 


W, + °° + dy Wel S —————; 
Ch i Im ml Snax lq" 


with integers q; (i =1,...,m) such that |q;| < B. The theory of diophan- 
tine approximations is concerned among other things with the study of 
this function A(B,c). For special numbers, nothing seems to be known 
about it except for quadratic irrationalities [1]. On the other hand, there 
is a theorem of Schmidt [2] stating that for almost all m-tuples, there is 
a number c, such that A(B, c) is asymptotic to c, log B. 

We have carried out computations giving this function for large values 
of B (up to 10* or 10°) and a few classical numbers like e, 7, log 2, y, 
e+, (/5 —1)/2. These computations in each case have a tendency to 
confirm the expected behavior, and we thought it would be useful to 
have them available in the literature. Most of them are for the usual 
case |qw + p| (i.e. m = 2). 

Of course, one knows the asymptotic theorem for ee — 1)/2 [1], but 
we have included it for comparison with the other numbers, about which 
nothing is known. 

One could ask for the behavior of another function t(B, c), equal to 


1 We are much indebted to Columbia University for use of the machine, IBM 794, Project 
No. UR7HBO1. 
Reprinted from J. reine angew. Math. 220 (1965), pp. 163-173. 
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the number of solutions of the inequality 
C 
lqgwy + <* +O, Wal S B 


with integers gq; such that |q;| < B. This is a more delicate function, and 
we merely observe that in all cases which we have computed below, its 
values are small, depending on c, ranging between 2 and 8, except in the 
case of 2, where there is a well-known disturbance corresponding to one 
unusually good approximation with gq = —113 and p = 355. 

The tables are easy to read. On the top of each table, we indicate the 
value for w,,w,... and c. In case we deal with q,w, + q2, we write it as 
qw + p. 

The columns from left to right give: 

The values for q,,q2,.-. (or q, Pp). 

In the third column, the maximum number N such that the absolute 
value of the sum is <c/N™™?. 

The absolute value of the sum itself. The symbols E-01 at the end of 
a number mean that this number should be multiplied by 107*. Simi- 
larly, E-02 means multiplication by 10~?, etc. 

Finally, in the last three columns, we give certain values for B, A(B, c), 
and A(B, c)/log B. Since c is fixed at the top of the table, we write A(B) 
instead of 4(B, c). Values of B were selected, for which a new solution for 
the inequality occurs. 

The quotient 4(B)/log B is the one which should be more or less 
constant. Its variations seem to be small, and also seem to follow some 
wave pattern. It is hard to tell whether there is anything significant to 
this. We have computed it by slide rule, and rounded it off to one 
decimal. 

If one plots the graph of A4(B) against B, one finds that it is a step 
function, which fits a curve c, log B rather well. 

For simplicity of programming, a few uninteresting solutions of the 
inequality in some cases have been omitted from the first four columns. 
These always occur at the beginning, with values of g,,q>,... equal to 
0, 1, 2, or 3. However, their existence was taken into account when com- 
puting the function A(B). 

We observe finally that we usually took a number c which is such that 
given an integer n> 0 there always exists a solution for an inequality 
(say) 


G 
lqw + p| <- 
n 


with |q| <n. This allows one to make a check on the computations to 
see that the required solution has been found by the machine. 
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w=e c=l+e 
q p N 1qw-+ pl ele 
log B 
—1 ? 5 7.1828182E-01 10 9 3.9 
-1 3 13 - -2.8171817E-01 20 11 3.7 
—2 5 8 4.3656365E-01 30 12 3.5 
—2 6 , 6 5.6343634E-01 40 13 35 
—3 8 24 1.5484548E-01 50 14 3.6 
—4 ’ ji 29 1.2687268E-01 70 US 3.5 
—7 19 132 2.7972799E-02 90 16 3.6 
—11 30 37 9.8899887E-02 110 17 3.6 
—14 38 66 5.5945598E-02 200 18 3.4 
—18 49 by? 7.0927087E-02 300 19 3.3 
—25 68 86 4.2954288E-02 390 20 3.4 
—32 87 248 1.4981489E-02 $80 72) 3:3 
—39 106 286 1.2991310E-02 1080 ap Bee 
—71 193 1868 1.9901794E-03 1270 P33} 3.2 
—110 299 337 1.1001130E-02 1460 24 3h) 
— 142 386 934° = 3.9803587E-03 2730 2S 3.2 
—213 579 622 5.9705381E-03 4180 26 Sh 
— 394 1071 1222 3.0404128E-03 $450 Pa 3.1 
— 465 1264 - 3540 1.0502334E-03 8170 28 3.1 
—536 1457 3955 9.3994594E-04 20510 29 2.9 
— 1001 2721 33714 1.1028750E-04 23230 30 3.0 
— 1537 4178 4481 8.2965844E-04 25950 31 3.1 
— 2002 5442 16857 2.2057501E-04 49180 32 3.0 
— 3003 8163 11238 3.3086251E-04 75120 33 3.0 
— 7543 20504 22141 1.6793342E-04 98350 34 3.0 
— 8544 23225 64502 5.7645921E-05 147520 35 3.0 
— 9545 25946 70633 5.2641582E-05 468490 36 2.8 
— 18089 49171 ' 743012 5.0043345E-06 517660 37 2.8 
— 27634: TSA 78054 4.7637251E-05 566830 38 2.9 
— 36178 98342 371506 1.0008669E-05 1084490 39 2.8 
— 54267 147513 247670 1.5013007E-05 ~ 1651310 40 2.8 
— 172346 468485 489086 7.6025026E-06 2168970 41 2.8 
— 190435 517656 1431122 2.5981571E-06 
— 208524 566827 1545336 2.4061301E-06 
— 398959 1084483 19357453 1.9208528E-07 
— 607483 1651310 1679450 =2.2139866E-06 
— 797918 2168966 9678726 3.8417056E-07 
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w=7, c=14+7 
A(B 
4 p N law + Pl BB 
—1 R) 29 =: 11.4159265E-01 10 10 4.3 
—1 4 4  8.5840733E-01 20 st 3.7 
—2 6 14 =2.8318531E-01 30 13 3.8 
—3 9 9 4.2477795E-01 50 14 3.6 
—6 19 27 ~—-1.5044408E-01 70 15 3a 
—7 22 467 = 8.8514247E-03 90 16 3.6 
—8 25 31 1.3274123E-01 340 17 3.0 
—14 44 233 1.7702850E-02 360 18 3.1 
—21 66 155 = 2.6554275E-02 380 19 Bes 
—28 88 116 3.5405699E-02 710 20 Sil 
— 106 333 469 8.8212804E-03 1070 2))\ Shi 
—113 355 137391 3.0144353E-05 1420 22 3.0 
—120 Si 466 8.8815690E-03 1780 23 Sal 
—226 710 68695  6.0288706E-05 2130 24 cel 
— 339 1065 45797  9.0433060E-05 2490 25 32 
—452 1420 34347 = 1.2057741 E-04 2840 26 33 
— 565 77/3} 27478 = 1.5072176E-04 3200 Zn Shs) 
—678 2130 22898  1.8086612E-04 3550 28 3.4 
—791 2485 19627 ~=—.2.1101047E-04 3910 29 S15) 
—904 2840 17173 =. 2.4115483E-04 4260 30 3.6 
—1017 3195 15265 2.7129917E-04 4620 31 Shi 
—1130 3550 13739  3.0144353E-04 4970 32 3.8 
— 1243 3905 12490 =—s- 3.3158 788E-04 5330 33 Shy 
— 1356 4260 11449 3.6173224E-04 5680 34 3.9 
— 1469 4615 10568 3.9187659E-04 6040 35 4.0 
— 1582 4970 9813  4.2202094E-04 6390 36 41 
— 1695 5325 9159 = 4.5216529E-04 6750 37 4.2 
— 1808 5680 8586 4.8230965E-04 104000 38 aa 
—1921 6035 8081 5.1245400E-04 104400 39 3.4 
— 2034 6390 7632 5.4259835E-04 208400 40 3.3 
—2147 6745 7231 = 5.7274271E-04 312700 41 aye 
— 33102 103993 216504  1.9129322E-05 521100 42 Br 
— 33215 104348 375994  =1.1015029E-05 625400 43 3.2 
— 66317 208341 510407 = 8.1142933E-06 833800 44 Sey 
— 99532 312689 1427780  2.9007206E-06 1146500 45 Shy 
— 165849 521030 794391 5.2135438E-06 1980200 46 3h 72 
— 199064 625378 713890 = 55.80 14411 E-06 2292900 47 a2 
— 265381 833719 1790708  2.3128232E-06 3126600 48 3:2 
— 364913 1146408 7044754 5.8789737E-07 
— 630294 1980127 2400864 = 1.7250422E-06 
— 729826 2292816 3522377 1.1757948E-06 
—995207 3126535 3642097 =: 1.1371449E-06 
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w=y, c=1+y 


q p N |qw + p| - B A(B) C2) 
log B 
—1 2 10 1.5443133E-01 10 9 3.9 
—2 3 5  2.6835300E-01 20 12 4.0 
—2 4 5 3.0886266E-01 30 13 3.8 
—3 = 13. 1.1392167E-01 50 14 3.6 
—4 7 38  4.0509654E-02 60 15 8h) 
—7 I 21.  =7.3412020E-02 80 16 Sh) 
—8 14 19  8.1019307E-02 100 il7/ i, 
—11 19 47  3.2902366E-02 130 18 Bui. 
—15 26 207 ~=—7.6072874E-03_ 150 19's 
—26 45 62 2.5295079E-02 250 20 3.6 
—30 52 103. = 1.5214575E-02 _ 280 21 ahi 
—41 71 89 —- 1.7687792E-02 400 22 BET 
— 56 97 156 1.0080504E-02 520 23 Si 
—71 123 637  =2.4732171E-03 790 24 3.6 
— 86 149 307 = 5.1340703E-03 1190 25 3:5 
— 142 246 318  4.9464341E-03 1580 26 aS 
—157 ~ Coie 592 2.6608532E-03_ - 4870 27 3, 
—228 395 8405 =: 1.8763610E-04 5260 28 3:3 
— 299 518- ~ 690 2.2855810E-03 5660 29 3.3 
— 456 790 4202 3.7527221E-04 10520 30 32 
— 684 1185 2801 5.6290831E-04 10920 31 3:3 
—912 1580 2101 . .7.5054441E-04 16170 a2 313) 
— 2807 4863 7117 =. 2.2158384E-04 
— 3035 5258 46460 3.3947739E-05 
— 3263 5653 10262 1.5368836E-04 
— 6070 10516 23230 6.7895478E-05 
— 6298 10911 13171 1.1974063E-04 
— 9333 16169 18383 8.5792886E-05 
w=et+n2, c=l+w 
A(B) 
q p N lqw + pl  &B A(B) jog B 
—1 5 7 8.5987448E-01 10 7 3.0 
—1 6 48 1.4012552E-01 20 8 all 
—2 12 24 2.8025103E-01 40 9 2.4 
—6 a0 43 1.5924689E-01 50 11 2.8 
_-7 41 358  1.9121374E-02 90 12 2:7 
—8 47 ' 56 1.2100414E-01 260 it} 2 
—14 82 179  3.8242748E-02 300 14 a 
—43 Psy 270 =2.5397272E-02 340 15 2.6 


— 50 Zs 1093 6.2758975E-03 630 16 203 
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A(B) 
q P N lqw + p| B A(B) jog B 
= 334 534 1.2845477E-02 920 17 2.5 
AGT 627 1044 6.5695792E-03 1840 18 2.4 
57 920 23358  2.9368167E-04 2760 + ~©=19 2.4 
= 314 1840 11679  5.8736333E-04 3680 20 2.4 
a7) 2760 7786  8.8104500E-04 4600 21 2.5 
—628 3680 5839 1.1747267E-03 19613 22 a2 
= 785 4600 4671  1.4684083E-03 20533 = s-23 23 
— 3347 19613 63176 — 1.0858254B-04 40146 24 2.3 
— 3504 20533 37060 —1.8509913E-04 
—6851 40146 89652 —7.6516590E-05 
w=loglog2+7, c=1+w 
A(B) 
q Pp N lqw + p| B A(B) ices 
=4 2 4 _7.7507973E-01 10 9 3.9 
acai 3 16  2.2492026E-01 208 11 3.7 
=9 5 6  5.5015946E-01 aor = 12 3.5 
=o 6 8  4.4984053E-01 - 40 13 3.5 
=3 8 11 3.2523920E-01 so) | 14 3.6 
=4 11 37 —-1.0031893E-01 Jon 15 3.5 
5 14 30 = 1.2460133E-01 90 16 3.6 
~9 25 155  2,4282403E-02 120 17 3.6 
=i 36 49  7.6036528E-02 140 18 3.6 
=18 50 717  4.8564805E-02 230 = «19 3.5 
05) 61 72  5.1754125E-02 340 = 20 3.4 
ll 86 137. 2.7471723E-02 700 ~=— («21 3.2 
—40 111 1183 3.1893203E-03 80 «22 3.3 
—49 136 178  2.1093082E-02 920 23 3.4 
~80 222 591 6.3786406E-03 1720 24 3.2 
Sly 333 394 9,5679609E-03 263025 a8 
—249 691 733  5.1464809E-03 4350 26 3.1 
— 289 802 1928  1.9571606E-03 6980 27 3.1 
—329 913 3063 =: 1.2321597E-03 11320 = -.28 3.0 
—618 1715 5206 —_7.2500097E-04 15660 29 3.0 
—947 2628 7443 5.0715869E-04 22630 ©. 30 3.0 
— 1565 4343 17329  2.1784228E-04 26980 31 3.1 
2512 6971 13048  2.8931642E-04 
—4077 11314 52817 7.1474140E-05 
— 5642 15657 25791  1.4636814E-04 
—8154 22628 26408  1.4294828B-04 
—9719 26971 50405 —-7.4893998E-05 
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w = log log 31, c=1+w 


q p N |qw.+ p| B A(B) MY 
log B 
a 2 2 7.6627796E-01 10° 12 5.2 
=) 2 4. 46744407E-01 20 «14 4.7 
=o 3 - 4 5.3255592E-01 30 «15 44 
3 3 , 3 7.0116611E-01 40 16 4.3 
=3 4 7 2.9883389E-01 - 60 17 42 
aay 5 34 = 6.5111854E-02 80 18 41 
= 6 13 1.6861018E-01 100 »3=—-: 119 4.1 
= 10 17 — 1.3022371E-01 140 20 41 
~9 11 21 =: 1.0349833E-01 7190221 4.0 
is 16 58 3.8386472E-02 30m 22° 4.0 
me Ph 83 2.6725382E-02  —«-:-330—St—s«é*F28 4.0 
—30 37 191  1.1661090E-02 420-24 4.0 
47 58 148  1.5064292E-02 650 =—-25 3.9 
—60 14, 95  2.3322180E-02 740 26 3.9 
= 95 656 3.4032016E-03 1070 27 3.9 
= 107 132 270 = 8.2578886E-03 1480 ~=—s- 28 a8 
SA 190 328  6.8064032E-03 - 1800 29 3.9 
—184 5B 460  4.8546870E-03 2540 30 3.8 
—261 322 - 1538  1.4514855B-03 3280 = 31 3.8 
— 338 417 1144  1.9517161E-03 4340 32 3.8 
— C0 644 _ 769 2.9029709E-03 5880 33 3.9 
—599 739 4465 5.0023061E-04 7620 34 3.8 
— 860 1061 2348  9.5125487E-04 10160 35 3.8 
1198. — 1478 2232  1.0004612E-03 
— 1459 1800 4952  4.5102425E-04 
— 2058 2539 45394  4.9206364E-05 
—2657 3278 4065  5.4943697E-04 
a5 4339 ° 5559 4.0181788E-04 
—4116 5078 22697  9.8412728E-05 
—6174 7617 15131  1.4761909E-04 
993) 10156 11348 — 1.9682546E-04 
wos a e=1, 0<qwt+p 
A(B) 
q P N Iqw + pl B A(B) jog B 
=) 0 1  6.1803398E-01 10 3 1.3 
=? 1 4 _2.3606797E-01 20 4 1.3 
—5 3 11 9.0169942E-02 40 5 1.3 
=i 8 29 + 3.4441853E-02 90 6 1.3 
—34 21 76 1.3155617E-02 240 7 1.3 
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AB 
q P N law +I BO 4B) “ J 
=o2 55 199 5.0249987E-03 610 8 i 
—233 144 521 1.9193787E-03 1600 9 12 
—610 377 1364 7.3313743E-04 -s- 4200_~—Ss 10 12 
— 1597 987 3571 = 2.8003358E-04 
—4181 2584 9349 =: 1.069633 1 E-04 
w, = log 3, w,=log2, ¢=log3 + log 2 
A(B) 
q Pp N lqw, + pw,| B A(B) ion 
—1 Z 6 2.8768207E-01 1 6 8.6 
—2 5 15: 1.1778303E-01 5 8 5.0 
=e) 2) 10 =: 1.6989903E-01 10 10 4.4 
—4 6 7  2.3556607E-01 20 13 44 
—5 8 34 5.2116001E-02 30s «14 4.1 
—7 11 27. 6.5667034E-02 40 15 4.1 
—10 16 17 1.0423200E-01 | 16 4.1 
—12 19 132 1.3551033E-02 70 17 4.0 
—17 P| 46  3.8564967E-02 90 18 4.1 
—24 38 66  2.7102067E-02 110 19 4.0 
27 46 71 = 2.5013934E-02 150 20 4.0 
—41 65 156 1.1462901E-02 0 621 4.1 
= 3S 84 858  2.0881324E-03 240 22 4.0 
— 65 103 114 1.5639166E-02 260 823 4.1 
—94 149 191 9.3747685E-03 320-24 4.2 
— 106 168 429  4.1762647E-03 a0 8625 4.2 
—147 233 245 —-7.2866362E-03 490 26 4.2 
—159 252 286 — 6.2643971E-03 | 4.3 
— 200 317 344 5.1985038E-03 1060 28 4.0 
—255 401 576 =: 3.1103715E-03 1540 29 4.0 
— 306 485 1752 —1.0222391E-03 2110 30 3.9 
359 569 1680 =: 1.0658932E-03 3170 31 EA) 
— 665 1054 41044 = 4.36541 10E-05 4220 a2 3.8 
—971 1539 1830 = 9.7858500E-04 5270 33 39 
— 1330 2108 20522 = 8.7308221E-05 6330 34 3.9 
Se 3162 13681 1.3096233E-04 23680 a5 3 
— 2660 4216 10261 = 1.7461644E-04 24730 36 3.6 
=P 5270 8208 2.1827055E-04 25790 37 oa 
— 3990 6324 6840 —_2.6192466E-04 50510 38 os) 
— 14936 23673 28970  6.1848688E-05 75240 39 35 
— 15601 24727 98477 = 1.8194580E-05  — 101020 40 3.5 
— 16266 25781 70376 = 2.5459532E-05 125750 41 35 


— 31867 50508 246630 7.2649564E-06 176260 42 a5 
— 47468 ta235 163936 —1.0929623E-05 302000 43 38 
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ee ee a en 
A(B) 
log B 


— 63734 101016 123315  1.4529913E-05 427740 44 3.4 
—19335 125743 488928 3.6646670E-06 478250 45 Sh) 
— UPt202 176251 497672 - 3.6002748E-06 603990 46 35 
— 190537 301994 27819424 6.4406776E-08 905990 47 3.4 
— 269872 427737 480483 —3.7290738E-06 1207976 48 3.4 
— 301739 478245 506736 =: 3.5358826E-06 - 1509970 49 3.4 
— 381074 603988 13909712  1.2881355E-07 
—S71611 905982 9277333  1.9313302E-07 
— 762148 1207976 6954856  2.5762711E-07 
— 952685 1509970 5566400 3.2188836E-07 


q Pp N-  ° |qw+pw,| ~ BAB) 


Wie log 33 W2 = log 5, W3 = log i a= Wight Ws 


3 
q1 q2 Ia COUN |sum| B A(B) a 
og B 
=) 0 1 2 2.5131443E-01 10 13 5.7 
2 at pe 4 — 8.5157807E-02 20 16 5.3 
at =3 3 4 — 8.9195578E-02 30 18 5.3 
= =a 19  4.1589966E-03 40 20 5.4 
=3 Aa 5 19  4.0377704E-03 50 23 5.9 
= =a ie 12 113 1.2122627E-04 60 25 6.1 
—6 —8 10 13. 8.0755408E-03 70 26 6.1 
S57 8 434 = 1171 1.1312253E-06 90 27 6.0 
= = 8 9 1.6272308E-02 150 28 5.6 
masks) =o 19 19 4.2802230E-03 ~—«- 180 29 5.6 
=i = 17 19 ©: 3.9165441E-03 220 30 5.6 
O02 =14 24. 79 2.4245255E-04 290 31 5.5 
— 33 =) ee 36 65 3.6367882E-04 360 32 5.4 
—44 = 28 48 56  48490510E-04 430 33 5.4 
=10 50 47 80 2.4150278E-04 700 34 5.2 
ot =57 59 113 1.2027651E-04 730 ° 35 5.3 
=) —64 71 1278  9,.4976534E-07 760 36 5.4 
=A3 5 83 112  1.2217604E-04 . 790 37 5.5 
—789 —72 505-2923 1.8146011E-07 830 38 57 
— 64 128 142 903 1.8995307E-06 860 39 5.8 
Bes i — 136 576 1420  7.6830519E-07 
—96 —192 213 737 -2.8492960E-06 
= 858 —200 647 950 —_‘:1.7180705E-06 
= 108 —256 — 284 639 . 3.7990613E-06 
— 160 = 320 355 571 4.7488267E-06 
—192 —9e4 426 $24 5.69859 19E-06 
3 = 0 4 — 7.6961040E-02 


34 —22 —1 46 7.2640789E-04 
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OOO 


A(B) 
q1 q2 q3 N |sum| B A(B) jog B 
5 —1 =2 13 8.1967671E-03 
693 — 120 —292 715 3.0307561E-06 
125 26 = 363 863 2.0809908E-06 
= 2 —1 ps 1.7435338E-01 
—1 3 —2 2 1.6211885E-01 
SL 43 S28) 65 3.6272906E-04 
= 2 = 11 1.2234538E-02 
=) 6 — 8 2.0431305E-02 
—12 36 See 56 4.8395534E-04 
es, 29 =i! 50 6.0518160E-04 
—174 149 =o) 189 4.3001952E-05 
w, = log3, w,=log5, w,;=log7, c=w,+w,+W, 
A(B 
Pe a er] BB) 
=) 0 1 4 2.5131443E-01 10 a2 neo 
= 1) 0 4 7 9.3354574E-02 20 41 13.7 
=e = 2 7 8.5157807E-02 30 46 iD 
—4 el 3 5 1.6615662E-01 40 a8 14.4 
—=( 7: 2 2 6.7294447E-01 50 57 14.6 
=" = 2 3 4.2566781E-01 60 62 15.1 
aa 4 3 3 4.2163004E-01 70 63 14.8 
=O} =) 4 5 1.7031562E-01 80 65 14.8 
—6 <2 : 7 8.0998811E-02 90 66 14.7 
Ste) = i) 19 1.2355764E-02 100 67 14.6 
= =3 3 7 8.9195578E-02 110 69 14.7 
—8 = us 33 4.1589966E-03 120 70 14.6 
a anh 4 4 2.4727665E-01 130 71 14.6 
—o 4 5 33 4.0377704E-03 150 i 14.4 
=) =) ie ii 8.1120037E-02 160 ie 14.4 
—0 6 3 7 7.2923270E-02 170 74 14.4 
= =6 14 23 8.3179933E-03 180 5 14.4 
=i 7) 12 195 1.2122627E-04 190 76 14.5 
=O Sy 10 24 8.0755408E-03 210 78 14.6 
= Jo) =, 434 2028 1.1312253E-06 220 8t 15.0 
= =" 8 16 1.6272308E-02 240 82 15.0 
—io =10 19 32 4.2802230E-03 270 83 15.0 
—14 — it 17 34 3.9165441E-03 280 84 14.9 
= —iZ 15 19 1.2113311E-02 290 85 15.0 
—4 —13 13 15 2.0310078E-02 310 86 15.0 
—22 —14 24 138 2.4245255E-04 340 87 14.9 
i —1 Ze 24 7.9543145E-03 360 90 153 
=o — 17 31 az 4.4014493E-03 390 91 15:3 
=20 —15 29 35 3.7953179E-03 450 92 152 
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qi . |Q q3 _N |sum| | B A(B) AB) 
log B 
—33 —21 36 113° — 3.6367882E-04 510 93 14.9 
—44 , —28 48 97 4.84905 10E-04 520 94 15.0 
22) 35 60 87 6.0613137E-04 540 95 15.1 
— 66 —42 Te, 79 7.2735764E-04 570 a7 153 
—270 —43 188 340 4.0152656E-05 600 98 153 
—10 —50 47 * 138 2.4150278E-04 630 99 15.4 
=r —57 = 196 1.2027651E-04 640 101 15.6 
—32 — 64 71 2213 9.4976534E-07 670 102 15.7 
—43 —71 83 195 1.2217604E-04 700 103 eed 
—789 —72 505 4236 1.814601 1E-07 730 104 - 158 
—54 +78 95 138 2.4340231E-04 760 105 15.8 
—65 —85 107 112 3.6462859E-04 790 106 15.9 
—31 —107 106 h13 3.6177929E-04 830 107 15.9 
— 302 —107 259 344 3.920289 1E-05 860 108 16.0 
—42 —114 118 - 139 2.4055302E-04 890 109 16.0 
—53 —121 130. 197 1.1932674E-04 920 110 16.2 
—64 —128 142 1565 1.8995307E-06 950 111 16.2 

—75 — 135 154 194 1.2312580E-04 

—821 — 136 576 2461 7.6830519E-07 

— 334 —171 _ 330: - 348 3.8253126E-05 

—96 —192 va) 1278 2.8492960E-06 

— 523 —200 - 647 1645 ~~ 1.7180705E-06 

—128 —256 . 284 1106 3.79906 13E-06 

— 885 — 264 718 1320 2.6678358E-06 

—160 - —320 © 355 989 —§ 4.7488267E-06 

—917 —328 789 1134 3.6176010E-06 

—192 — 384 426 903 5.69859 19E-06 

—949 —392 - 860 1009 4.5673665E-06 

—224 —448 497 836  6.6483572E-06 

—256 —512 568. 782 7.5981227E-06 

—288 —576-- 639 dod 8.5478879E-06 

2 —1 0 Z 5.8778666E-01 

2 —2 0 2 1.0216512E-00 

3 —1 —1 4 2.5951119E-01 

12 -—7 —1 12 2.8628072E-02 

34 =—22 =1 80 7.2640789E-04 

5 —1 —2 23 8.1967671E-03 

8 2 —2 8 6.8764273E-02 

10 —2 —4 16 1.6393534E-02 

37 —18 —6 37 3.3113625E-03 

565 — 376 =o 825 6.8298 175E-06 

45 —15 13 74 8.47634 16E-04 

56 —8 —25 69 9.6886044E-04 

206 - —85 —46 aaZ 4.2052187E-05 

217 —78 —58 242 7.9174088E-05 


597 —312 —79 889 5.8800520E-06 
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4(B) 
91 q2 q3 N |sum| B A(B) log B 
238 —21 —117 336 4.1102422E-05 
629 — 248 —150 971 4.9302868E-06 
661 — 184 —221 1081 3.9805214E-06 
693 —120 — 292 1239 3.0307561E-06 
qa — 56 — 363 1495 2.0809908E-06 
—1 2 —1 5 1.7435338E-01 
-1 .3 —2 5 1.6211885E-01 
-1 43 — 35 113 3.6272906E-04 
—2 a —1 2 9.2425889E-01 
—2 5 —3 19 1.2234538E-02 
—3 3 —1 3 4.1343327E-01 
—4 4 —1 6 9.7392346E-02 
—4 10 —6 13 2.4469075E-02 
—5 7 —3 8 6.4726502E-02 
—7 6 —1 iS) 2.0431305E-02 
—9 11 —4 tht 3.2665842E-02 
—12 36 —23 98 4.8395534E-04 
—15 32 —18 36 3.5538151E-03 
—23 29 —1li 87 6.0518160E-04 
— 26 WS —6 36. 3.4325888E-03 
—31 26 —4 31 4.7641783E-03 
—35 65 — 34 65 1.0891369E-03 
—46 58 —22 62 1.2103632E-03 
—57 Sil —10 59 1.3315895E-03 
—110 DE — 167 321 4.4901483E-05 
—142 213 — 96 325 4.3951718E-05 
—153 206 — 84 245 7.72T4557E-05 
—163 156 —37 168 1.6422823E-04 
—174 149 —25 328 4.3001952E-05 
—185 142 —13 243 7.8224322E-05 
—327 355 — 109 368 3.4272604E-05 
—359 291 — 38 363 3.5222370E-05 
— 437 632 —276 661 1.0628878E-05 
— 469 568 —205 693 9.6791134E-06 
—501 504 —134 730 8.7293481E-06 
—533 ~ 440 —63 773 7.7795827E-06 
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APPENDIX B 


Continued Fractions for 
some Algebraic Numbers 


By S. LANG and H. Trotter at Princeton 


The following tables at the end of the paper contain the continued 
fractions and some related information for a few algebraic numbers, viz: 


m2, Vs iE oS. vt Baie 2 ies 


where 
2n . Cee 
a, = 2 cos-, is a root of x” + x* — 2x — 1, 
a, is a root of x° —x —1, 
a, = 3/2 + ,/3 is a root of x® — 9x* — 4x3 + 27x? — 36x — 23. 


The first numbers were chosen as cubic irrationalities, «, as being totally 
real, and we picked «,,a, at random to be non-cubic. The original 
motivation was to see if computations for algebraic numbers would be in 
line with some conjectures made in [3] and [4], and to observe whatever 
else might come up. 


APP. B, 1. TABLE | 


In each case, this table is to be read horizontally, and gives the first 
thousand terms in the continued fraction. For instance, for eo the con- 
tinued fraction begins 


(iver), 1, ty 40 s,s. - 


Reprinted from J. reine angew. Math. 255 (1972), pp. 112-134. 
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The last five terms among these first thousand are [..., 2,1, 1, 1,1]. Thus 
the continued fraction [a,,a,...] is laid out in rows of twenty integers, 
and there are fifty such rows. The position of any a, is therefore easily 
determined. 

For clarity, and reasons of space distribution, the table is filled only 
with the 2-digit integers which arose. Whenever a larger integer oc- 
curred, it is indicated by the symbols *a, *b, *c,..., and these are then 
listed separately at the bottom. For instance, for 2 we have 


*q = 534, *h = 121, seis *k = 4941, 


APP. B, 2. TABLE Il 


This table gives the frequency count for the digits appearing in the 
continued fraction. According to a theorem of Kuzmin [2], for almost 
all numbers «, the probability that the n-th integer a, in the continued 
fraction for a is equal to a positive integer k is given by 


(k + 1)? 
a 1) 

For k = 1, and almost all numbers, this means that the probability for 
a, = 1 is approximately 0.41. In each case of the computations, we find 
that the numbers behave very closely to this generic expectation, thus 
confirming that in most respects, the continued fractions for algebraic 
numbers of degree >2 should be essentially like those of almost all 
numbers. For instance, among the first thousand a, for ee, 0/4, and a, 
we find that 1 occurs, respectively, 422, 412, and 418 times. The biggest 
divergence from the Kuzmin number is for 5; and even then 433 is still 
in line with the expected asymptotic estimate. The Kuzmin probability 
that 2 occurs is approximately 0.17, and again in this case we find that 
2 occurs 165 times (resp. 159 times, resp. 168 times) for 2, 8, 
a%,, respectively. One should also note that this type of regularity is 
exhibited throughout, among the first thousand terms. 


APP. B, 3. THE OCCURRENCE OF LARGE INTEGERS 


Aside from what appear to be routine low numbers, there occur larger 
numbers which seem to be of two kinds: Some which are only some- 
what larger than the ordinary ones, and some which appear to be excep- 
tionally large. For instance, for 3/2, we meet 


a572 = 7451 and 4620 = 4941. 
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For 3/5, we meet 
Aig = 3052, 4691 = 13977, - agi3 > 49968. 


The occurrence of large numbers in the continued fractions of certain 
cubic irrationalities has already been observed by Brillhart, see Church- 
house and Muir [1]. A _ theoretical explanation for some of them 
has been proposed by Stark [6], but the following problems remain: 
Determine whether there is a basic theoretical distinction between what 
seem to be only medium large numbers, and very large ones. More 
importantly, determine whether exceptionally large integers will continue 
to occur throughout the continued fraction, or whether they will stop 
from occurring. The explanation given by Stark depends on some class 
numbers being equal to 1, and thus would account for only a finite 
number of them. In general, the appearance of such large integers may 
depend on the arithmetic properties of the field obtained from the square 
root of the discriminant, e.g. its class number. The tables seem to indi- 
cate that they stop. 

To discuss the statistical significance of exceptionally large values of a, 
occurring near the beginning of the sequence of partial quotients, we 
need an estimate of the probability qy , that the first N partial quotients 
of a “random” number are all less than a given integer K. It is perhaps 
most natural to consider a random number as distributed uniformly on 
(0, 1), but in this context the distribution given by 


Pr{X Sc}=log,(i+c) for ce(0, 1), 


is more appropriate, because if X has this distribution, then the distribu- 
tion of the partial quotient a,(X) is independent of n. In fact, 


Pr{a,(X) < K} = yx, 
where 


1 
w= Prix <K} =1— loea(1 a5 z) 
which is the Kuzmin theorem already alluded to. 
To see that this is so, observe that a,(X) = [X,'], where X, = X, and 


for n > 0, 
Xn4i = Xe aan ee: 


As usual, [x] is the largest integer < x. It is then an exercise in calculus 
to show that if f, is a density function for X,, so that 


Pr{X, Sc} = [, f00 dx, 
0 
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then f,., = Tf, is a density function for X,,,;, where T is the linear 
operator on L'(0, 1) defined by 


2) 


(T(x) = Y w+ 7f((x +). 


k=1 
It follows that f,,, = T"f,. It is easy to verify that the function 


1 
(log 2)(1 + x) 


is a density function and is invariant under T, so that if X has this 
density, in which case, by integration, 


Pr{X Sc} =log,(1 + 0), 


all the X, have the same distribution. In fact, Kuzmin’s theorem states 
that if f is any smooth probability density, then 


. : _ 1 
hoy SO) (eee! 


It follows that as n— oo the distribution of a,(X) tends to the one given 
above if X has any smooth distribution. For a discussion of all these 
ideas and a proof of Kuzmin’s theorem, see [2]. 

If the random variables a,(X) were independent, then we would have 


Qn,K = YR: 


This is not strictly correct, but can be expected to give a good approxi- 
mation for large values of N and K. A combination of theoretical and 
numerical analysis indicated strongly that the relative error is bounded 
by ANK~?, with 2 <1, and we are confident that the approximation is 
entirely adequate for our purposes. 

The short Table A at the end of this section shows, for each of the 
numbers investigated, the maximum value A of the first 1000 partial 
quotients, and the value M for which ay = A. The third column gives 


P1000,4 = | — 41000,4> 


the probability that a “random” number would have a value as large as 
A among its first 1000 partial quotients. The smaller the value of p, the 
stronger the evidence that the number is unusual. The fourth column 
gives 


Pu,a=1—4m,a> 
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the probability of getting a value as large as A among the first M 
quotients. The value is of course smaller than that of p, and its statistical 
meaning less clear since M is taken a posteriori to make the probability 
small. 

The table shows that if one goes by the maximum quotient found, 
only gS appears highly unusual, although one might question oA and 
a. If one also takes into account the second largest quotient, then v5 
with dg9, = 13977 appears even more unusual, and 2 with dg.9 = 4941 
perhaps comes to be of interest also. 


Table A 
Number A M_  P=1-—41000,.4 Pu =1—4m.a 
6) 7451 572 0.18 0.10 
A/S 3502 916 0.34 0.31 
3/4 14902 579 0.09 0.05 
J/5 49968) 813 0.03 0.02 
/7 689s G1 0.88 0.72 
oy 904 830 0.80 0.73 
Oy 11644 588 0.12 0.07 
3 1446 54 0.63 0.05 


APP. B, 4. TABLE Iil 


In each case, this table begins with the columns labeled n, a,, and q,. 
The n indicates n-th position in the continued fraction. The a, means the 
n-th partial quotient. The q, means the denominator in the approximating 
fraction p,/q, (classical notation). For instance, in the case of ia we 


have 
36 => 534, 136 = 3.06 x 10+% 


A4g2 = 121, qa2 = 8.95 x O72, 


and so forth. In machine language, E 19 means multiplication by 101°, 
and E 486 means multiplication by 10*°° (the last line in Table III). 

Table III includes these data for all n among the first thousand such 
that a, => 50. We picked 50 as a cutting point after looking at prelimi- 
nary computations, because it included all the numbers a, which could 
be labeled as somewhat large, and at the same time provided only a 
rather small table. 

The last column r, in Table III gives (up to three decimals) the quo- 
tient 

Qn 
Gn—1 108 Gy-1 
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for those values of n when a, = 50. The reason for this quotient to be 
interesting are as follows. According to a theorem of Roth, if « is alge- 
braic, there is only a finite number of integers q > 0, p such that 


1 
te Pl ara 


It was suggested in [3] and [4] that this theorem should be improvable 
by an inequality 


it 
ee Pl aie 


where f is a function close to the logarithm, for instance (log q)'*‘, or 
perhaps even log q itself, up to a constant factor of course. Such a 
function is called a type in [4]. If some a, is small, then q,/q,-, being 
approximately equal to a, shows that the quotient 


Qn 
Gn-1 108 n—1 


is approximately like 1/log qg,_,. Thus to investigate the possibility of a 
type f, we look at those n for which a, is comparatively large. Again n 
such that a, 250 seemed to give the most information for the least 
amount of space used. It is even unsolved for any algebraic number of 
degree >2 whether it is of bounded type, but the tables seem to fall 
fairly well in line with expectations, e.g. differing from the log by a 
function with a lower order of magnitude (above or below). 

We have also programmed the same data for the first 3000 terms of 
the continued fractions of the cubic numbers listed. In every case, excep- 
tionally large integers did not seem to recur, and generally speaking, the 
ratio r, seems to decrease. We thought it pointless to reproduce these 
more extensive tables in full, but we give in Table B the portion of Table 
III for n > 1000 when r, > 1, rounding off r, to one decimal. 

The tables therefore suggest that the type may in fact not be bigger 
than a constant times the logarithm, and may even be of an order of 
magnitude smaller than the logarithm. Following certain asymptotic esti- 
mates of Adams, who looked at the continued fraction of e, it was shown 
(cf. [4]) that the type of e is asymptotic to log g/log log g. Thus one is 
beginning to be accustomed to such small types. Note that for a func- 
tion essentially not bigger than the log, the series 


1 
L@ 


diverges, so that these cases go very slightly against the Khintchine 
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convergence principle: If y is such that )’ w(q) diverges (resp. converges), 
then for almost all numbers, the inequality 


|qa — p| < W(q) 


has infinitely many solutions (resp. only a finite number). However, this 
statistical result is delicate to use for specific numbers with a type in the 
range of the log, because one also knows that if w is a number such that 
for every function y (decreasing) having convergent sum, the above in- 
equality has only a finite number of solutions, then « must be of bounded 
type. For all this, cf. [4]. 


Table B 
n ‘ a, Qn Vn 


2 1191” 12737 7.74E1010 5.5 
2248 «=: 2897: 2.97E1136 = 11 

33 1988 2967 347E1024 13 
2407 «(9559 1.25H 1242 3.3 

3/4 1974 6368 488E£1010 2.7 
_ 2248 «= 4157 6.92 E1146 1.6 
aS 1196 18905 147E600 13.8 


ye COs 1102 = 1374S 6.84 E576 1.0 
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The computations were done by the following algorithm, which uses 
integer arithmetic only, and thus involves no rounding error. 

Given a polynomial p,(x), of degree d, with positive leading coefficient 
and a unique positive root y, which is simple, irrational, and greater than 
1, we construct a polynomial P,,,(x) as follows. Let a,=[y,] be the 
greatest integer such that P,(a,) <0. Define 


Q,(x) = P(x +a,) and — Pay (x) = —x7Q,(x~"). 


Then Q,(x) has exactly one root between 0 and 1, and since the roots of 
P.4,() are the reciprocals of the roots of Q,(x), we see that P,.,(x) has 
a unique positive root y,4,. Obviously y,4, is also a simple root, 
irrational, and greater than 1. Note that the constant term of Q,(x) is 
P,(a,) <0, so that P,,,(x) also has a positive leading coefficient. Thus 
P,,(x) has the properties assumed for P,(x), and starting from any P, (x) 
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with these properties, we can define an infinite sequence P,(x), P,(x),... 


with associated positive roots y,, y2,..-- 


We have 


Viet = (Vn aa Ga 


and 


a, = [yn] 


If P,(x) has 


integer coefficients, then so has every P,(x), and the calculation involves 


only addition and multiplication of integers. 


is the sequence of 


partial quotients in the continued fraction expansion of y,. 


This is precisely equivalent to saying that a,,4q),... 


Table I 
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Table I (continued) 


Gees IMO ee 1 203 4 Se 6 7185910 


2 13) 18, 5 


1 


1 


4 
22134 


600 
620 
640 
660 


2 


1 


1 


1 


5 


a2 (SS 49 66 2 
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a= 534 
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Table II 
Re 
Frequency Counts 

1 422 Ze 4 56 2 

2 165 23 3 67 1 

3 91 25 3 69 Zz 

4 66 26 2 71 1 

5 40 ef 4 de 1 

6 37 28 1 77 1 

gi 20 29 2 84 1 

8 16 30 1 89 1 

9 15 31 1 93 1 
10 15 33 1 108 1 
11 8 34 2 113 1 
12 9 36 1 121 1 
13 8 38 2 151 1 
14 a 39 2 186 Ds 
15 5 41 2 220 1 
16 3 44 1 255 1 
17 2 47 1 372 1 
18 3 48 it 534 1 
19 3 49 2 4941 1 
20 4 50 1 7451 1 
21 3 53 1 
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Table II 
n a, “- dy ' 
36 534 3.06E 19 13.844 
42 121 8.95 E 22 2.530 
73 . 69 448E 39 0.799 
89 89 3.61 E 48 0.835 
92 186 1.56E 52 1.618 
fit D2 S72 6.69 E 65 2.560 
269 72 1.90 E 146 0.219 
376 186 4.78 E 194 0.421 
421 220 © 4.51 E216 0.447 
468 84 1.19 E 239 0.154 
494 56 LOMB 253 0.098 
511 255 3.59 E 260 0.430 
528 71 3.38 E 268 0.117 
540 67 1.32 E 276 0.106 
Si 7451 8.64 E 297 11.005 
587 113 4.07 E 308 0.160 
611 69 7,64 E 320 0.095 
618 Sy ES 26 0.202 
620 4941 2.97 E 330 6.568 
654 33 1.50 E 347 0.068 
761 56 Sle} 13, hs; 0.063 
786 50 7.54 E 406 0.054 
798 77 4.14 E 414 0.082 
937 108 . 3.81 E 475 0.099 
956 93 1.37 E 486 0.084 
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Table I 
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Table I (continued) 


ae 
De see 910) ol 23 4 5 6. 7 8 G10 
S00ne ol) 42139 1°19 fi 5 1 #38 122128 sé5 3 1«1 
Pomel > Jeet 2 5 2 if 5 3 12 t 6-4 41 4 
mers 15 ol iets G3 i OS 1 OL 1 SA ES 
coed) 1 20 1 1 “to 2 2 21 112 262. 1 “hl Steed 
SSIS cr Ry ns a es Sy 
Som 1 610 8 4 1371 25 «131620 5 % 1 1 «7 8 
mms, fo 2 74 83) 4 2 12 2 11 23 2 2 
coats S ft £2 47 £ 11 ~«2 1 2a 2 2 Sho tat 
oemeeice 2 28 2 3 $0313 4 1 = 1 1 6392 £ «13 8 £2 1 
cegmee? 20501 6 PLS 1 %e° «1 «2 4-4 3 15 7 2 30 2 
a=139 b=249 c=612 d=220 e=123 f=131 
g=196 h=729 i=164 j= 396 k=343 m= 137 
n=139 p=268 q=247 r=1232 s=116 t = 3502 
u = 164 
Table IT 
a3 
Frequency Counts 
1 425 13 4 bi 46 1 84 1 247 1 
2 159 14 2 26 «3 os 1 s9 (1 249 1 
= 97 1s 8 pga 1 54 3 oF 1 268 1 
4 64 16 5 28 4 56 1 95 1 343 1 
5 37 7 2 31 2 a 1146 1 396 1 
6 32 18 3 Bay 1 60 3 123 1 612 1 
7 14 19 3 34.3 66 1 a | 729 1 
8 18 20 5 ar 62 67 1 i379) 1 123241 
9 10 a 4 38 69 1 139 2 3502 1 
10 13 22. 1 39 3 Tiel 164 2 
11 12 23. «4 Ate 2 I A 196 1 
12 9 24 4 a. 1 Pa | 220 1 
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Table III 
3 

n a, . Qn T, 

64 139 6.85E 30 2.118 

81 52 9.42E 38 0.614 
109 249 3.42E 54 2.077 
119 612 8.70E 60 4.575 
146 60 7.46E 73 0.363 
163 58 1.94E 84 0.307 
184 82 218 E 96 0.379 
194 54 3.33 E 102 0.233 
207 220 3.57 E 110 0.885 
PAU 89 9.94 E 113 0.347 
243 67 3.00 E 128 0.231 
253 60 4.08 E 133 0.200 
308 123 1.53 E 157 0.345 
310 84 1.31 E 159 0.235 
354 95 ~ 91996 E77 0.236 
371 131 9.09 E 185 0.310 
378 : 196 7.22 E 191 0.450 
386 729 1.24 E 197 1.631 
405 164 4.22 E 207 0.347 
477 66 2.34 E 242 0.121 
535 58 2.72 E 274 0.093 
580 396 9.29 E 296 0.585 
621 54 5.56 E 320 0.074 
631 92 2.84 E 326 0.124 
708 71 2.10 E 363 0.085 
719 343 5.43 E 369 0.406 
722 137 1.49 E 372 0.161 
731 139 4.61 E 376 0.162 
736 268 3.91 E 381 0.307 
751 247 3.47 E 389 0.278 
761 69 2.10 E 395 0.077 
766 1232 6.02 E 399 1.349 
816 56 2.01 E 425 0.057 
857 54 6.13 E 443 0.054 
896 - 116 1.77 E 462 0.110 
899 60 3.23 E 464 0.057 
916 3502 3.38 E 477 3.209 
974 ne 1.22 E 507 0.062 


990 164 LSE sts 0.139 
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Table I 


Coie) 9 10 91-2354 5 6 7 $9.10 


23245 °5 


1 


29 6 4 


=i 


40 


1,1 


1 45 


ieee Ss «62 10° 6 1 
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80 


5 
i172 


1 
1 


| 2 RS A Se Rs 
6 23 99 


ee 1th ok 4A 
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120 


4 


422 2 11 


P22 2 


160 


180 


11 3 
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200 


2 


leecorom Ooo! 1 18 1 “Siveaees 8 1 °3 eee 


13 


240 
260 


1 


Bee S 3 


1 


3 


i 


ise? 2 18 


2a 6 
453° 3 


: 
1 


1 
2 


280 
300 
320 
340 
360 


O15. (See 
2° Seon 0 
6 5 Gio 


oe eG" ..2 3 


1 


1 3? 


See 23. 


380 


ad 


1 
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Beet Seu f 


2 4 
2 
1 
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4 10 


PIGS bt I? 4 


1 


1 


400 
420 
440 
460 
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3 7 Taesels 


L167. 5 
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ieg2e 27 
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OM 


3 


22 3. oh 


1 


3 


14 


Omer! 
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500 
520 
540 
560 
580 
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1 
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PGS: .2.55 
58 2 6 4 34 
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600 


3 4 2 
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Seo 2 


640 
660 
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108 CONTINUED FRACTIONS FOR SOME ALGEBRAIC NUMBERS [APP. B, 5] 


Table I (continued) 


4 
1 232.5 faa. & Jaw le 3 UW 8 P67 Sop 
800 1. tae 2 or “Aes so Se 2 £ 4 PS ies ae 
820 $ 1.2310 2-8 2510 i492 5 3 te 2a 
840 i ot 1 1 eee ede s 3 2.2 § 2 bt Ae 
860 2 3 2 @ eile 4g a ot 2 eee 
880 1 is ee ee ee 1 2@ 146 3. 20 159i >2 8 
900 Talo 3) 74) ORB IEe I) 52a 1d of @ ae & 3e58.2 
920 1111 4a 3 154 3° 2, 2 2 to. 4 4 oe Wy al 
940 3 1 ad 44°23 Bee 1 88 ty Bee Ss 1 l 235 
960 3.1 1°96 19eia 1 Tes 522082 14i11f 1 ae 
980 21264 339 11 «21 ive 1 2 4 3 3 Gee 
a = 266 b = 745 oie d—at0 eS sili f = 144 
g = 14902 h = 139 i = 303 j = 2470 
Table II 
a 
Frequency Counts 
1 412 2D 1 54 if 
2 164 23 3 55 1 
3 100 24 5 56 2 
4 69 2S 2D 58 1 
5 39 26 2 60 1 
6 52 27 3 78 2 
7 19 28 1 79 1 
8 20 29 1 93 1 
9 14 30 1 99 1 
10 12 31 2 110 1 
ital 11 33 Z 139 1 
12 is 34 2 144 1 
13 6 35 1 266 1 
14 4 37 1 303 1 
15 4 38 1 Sp 1 
16 3 39 1 Sul 1 
17 4 41 3 745 1 
18 2 . 44 3 2470 1 
19 1 45 1 14902 1 
20 8 46 2 
21 1 50 1 
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Table III 
4 
n a, z= “an Th 
41 266 193E 19 6.868 
47 60 a7, 22 1.248 
101 93 985E S51 0.808 
124 745 8.43E 65 5.136 
145 78 ° 2.00E 76 0.460 
158 99 5.29E 86 0.508 
187 56 1.53 E 100 0.249 
318 78 2.14 E 167 0.205 
379 372 6.02 E 194 0.842 
404 79 1.26 E 207 0.168 
420 110 2.84 E 216 0.223 
508 511 4.53 E 260 0.861 
527 144 4.29 E 268 0.236 
565 5D 1.14 E 289 0.084 
579 14902 1.09 E 298 22.025 
581 58 6.43 E 299 0.086 
596 56 2.54 E 308 0.079 
618 139 9.69 E 320 0.191 
625 5 Sie (E26 0.407 
627 2470 1.87 E 330 3.283 
702 50 2.38 E 367 0.060 
926 54 2.40 E 475 0.050 
Table I 
¥S5 
fe? seed | Oar 8.9 10 feces 4 See 7 8 0 
0 Site 2 4. aa ft 5 1 1-410) 17 fle 1- 1 Xa, 1 
20 plete tt 1s Se bd: 2 PAs aie t ees: 2. 40 
40 Anizowle os 2° 7 tty 3: 33 Cee reeset 13 2.2 
60 teise 8§ 10 48 1.2 1 1 sea oO 1 132 4 
een 49° 310 1 8 1 1 1 1 4 imGume 2 2 2 3 “Seal 
100 Seo cemleol.~t 10° 2 °°§ 3.4 Ayes 1 CL OL. Tee 
120 2 wee? | 1 G1 3 1 Dias) /5 a) Gl Gee Ga ee 2) 
140 PIE 08 i) 2 a Iles ee a | BUI 45 3: 73) 61 GS 
160 DAMEN iis eae eas | se 2 a leat 1 See 66 oe) 
180 pao 2- FO a tle aan 2 S26, bee 2 bes 1 ooet 
Omid 15720 ft +242 1 1 2 So lesemes) | (22 SL 
promeieaemoe tl 52> fe2 1 2°92 5 413 2 3 1 1 1 9 
240 Pee eto) 6h 2 t2lhUe8 he dl dE wi 
260 ine te oe ies tf 917 16 ft 1 LS 2 %¢ 2 
280 Set edeete tt 25 2 oS 6 Bol <2 9273 4.1.49 4 
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Table I (continued) 


1 2535:4 5506. Ted OR 


6 7 8 9 10 


pee er, ee 


1 


300 
320 


1 38 


2 


1d 1. 3°*er 23 Ob Sista 
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360 
380 
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420 
440 


20 5 4 2 


1 


Deel t2 2 See Pt 


mw 2102 9 


480 
500 
520 
540 
560 
580 


125° Sipeeo 


oe oe tOn 2 


3 


1 


PI: 


ee 
ys 


6 6 


600 
620 


2 15 552562 


(Nes a | 
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1 


2 2 7a 


1 


128 20 “3 


3 


1 


640 
660 
680 
700 
720 
740 


1 t0P 


L 3 7 a2 


25 2 6 4 


5) oes 2 


780 


820 
840 
860 
880 


900 


1 


1 


9 247 


jh Sh tb a 327 +8? 
*n 13 


1 


1 


a2 fe 


$2 
6 2 


1 


920 
940 


4 


3 B23), 4 22 


980 


1051 


c = 854 
fas 16S 


= 474 
326 


b 
h 


@— 5052 
g = 182 
oo 


49968 


451 


k= 


j = 13977 


p= 121 


[APP. B, 5] 


COMPUTATIONAL METHOD Ltt 
Table II 
fs 
Frequency Counts 

1 433 ily 3 41 1 US ii 
2, 180 18 ile = 42 2 81 1 
3 95 19 1 - 45 2 IDA 1 
4 50 20 3 47 1 131 1 
5 46 Py| zt 48 2 135 1 
6 DS 23 4 49 1 170 1 
7 19 25 3 50 1 182 1 
8 19 2b ele 51 D 326 1 
9 16 29 2 - 52 1 45) 
10 12 Bil 1 53 2 474 1 
11 4 33 2 55 1 739 1 
12 8 34 1 59 1 854 1 
13 13 35) 22 60 1 1051 1 
14 2 ° BT 3) 61 2 3052 1 
iS 4 38 1 63 1 13977 1 
16 2) 40 2 72 1 49968 1 


112 CONTINUED FRACTIONS FOR SOME ALGEBRAIC NUMBERS [APP. B, 5] 


Table III 
Ss 

n a, Qn tT, 

19 3052 428E 11 162.730 

74 474 1.23 E 40 Syollé 

92 60 9.82E 49 0.548 
104 61 2.43E 56 0.491 
215 81 1.88 E 109 0.331 
226 52 2.20 E 116 0.200 
2A] 51 3.15 E 140 0.162 
279 854 5.44 E 143 2.636 
Bi 131 1.57 E 183 0.315 
31) 170 1.07 E 186 0.402 
378 51 1.34 E 189 0.120 
388 59 4.04 E 195 0.134 
406 1051 1.00 E 205 2.260 
415 61 2.21 E 210 0.127 
457 182 2.35 E 230 0.348 
464 50 6.48 E 233 0.095 
487 326 5.54 E 246 0.581 
495 135 8.42 E 251 0.236 
520 SS 9.33 E 265 0.091 
577 72 2.06 E 291 0.108 
691 13977 2.20 E 345 17.792 
762 75 1.46 E 379 0.087 
765 63 7.47 E 381 0.073 
809 451 4.93 E 401 0.491 
813 49968 1.73 E 407 53.913 
847 53 3.15 E 423 0.055 
871 739 8.92 E 438 0.737 


B27 121 7.42 E 462 0.114 
970 aS 9.37 E 484 0.048 
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[APP. B, 5] 


Table I 


ie 3 4 5686" 7) 8 9-10 


ca SS? «8 9 10 
Pit0 2 16592 
3 1 2 27 4 
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3 2 624 
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Seemed 22 See eee, 2 
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See 
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eed 3 23 I 
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Bees 2 63 (CS 
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580 
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640 
660 
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700 
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1 69 


1. Naas | 
Peis 2 2 
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7 


We) 
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2 4) 2 13-11 


2 
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*n 11 


Soe 


3 


114 CONTINUED FRACTIONS FOR SOME ALGEBRAIC NUMBERS [APP. B, 5] 


Table I (continued) 


ai 
12)).3. 4 SeerGaey 879510 1°2-3) 45.6 Tes 
B00 GS 1 Ct CS PS 24 1 ees PO 4 ee eee 
820 52 Toa 3 138 lotol 2 oe? 12 ee 
$400924 2 1 44 13° °2 £ 3 Boi 113 Pee 
860 A 9 2 42) 34 403 7 211 2 1 4 8333570 
SeOn olga. 1 Zee pads? 230k 15642 1 10) “i See see 
900 822 43 1 444 6 1 7°22 9 218 See 
920 24 2a 2 2 fo2 5 1 2 52 35s 
940 5 eel. l- Goma 3 24 16 ft 2 tela 
960 pS Ee Soi 2.4 2 *%q 2553 
OsGmecs 4 1-1) 2 2.12 3, ial ied. 4122 1) 9 ere 
a= 282 b = 104 c= 2/7) d = 429 e = 303 f = 341 
g = 110 h = 197 i= 118 j= 133 k = 689 m= 115 
n= 111 p = 202 q = 628 
Table IT 
ei 
Frequency Counts 
1 409 18 5 43 1 104 1 
2 161 19 7 44 zZ 110 1 
3 88 20 1 47 1 111 1 
4 55 21 2 51 1 115 1 
5 34 22 Zz a2 2 118 1 
6 40 23 2 53 1 133 1 
a 29 24 2 56 Z 197 i 
8 24 25 3 60 2 202 1 
9 13 27 a 63 1 227 1 
10 12 30 2 69 1 282 1 
11 sl 31 3 da 2 303 1 
12 ao 32 1 73 1 341 1 
13 8 33 1 78 1 429 1 
14 : 8 34 2 84 1 628 1 
15 11 35 3 87 1 689 1 
16 3 36 1 98 1 
17 4 41 1 99 1 


[APP. B, 5] COMPUTATIONAL METHOD He 
Table II 

a a, —~- Gy is 

29 282 6.14E 15 9.209 

89 104 2.74E 44 1.66 
116 Seel TI9E 59 2.096 
118 429 3.10E 62 3.120 
120 Te 2.27E 64 0.507 
262 56 « WHAs 18 135) 0.185 
310 303 4.11 E 160 0.834 
322 56 7.57 E 167 0.147 
370 341 2.07 E 195 0.770 
435 110 2.82 E 226 0.214 
444 197 . 146 E 232 0.374 
446 60 6.16 E 234 0.112 
458 98 4.21 E 242 0.178 
466 118 2.06 E 246 0.211 
542 60 2.86 E 283 0.093 
Soy jes, 1.21 E 293 0.200 
567 Doz 4.71 E 298 0.077 
611 689 P2008 321 0.940 
653 he) 6.74 E 339 0.127 
658 51 1.79 E 344 0.065 
679 69 1.36 E 355 0.086 
705 78 9.95 E 369 0.093 
Mo! 115 150 E373 0.135 
786 111 4.29 E 411 0.118 
801 63 9.92 E 421 0.065 
811 53 5.06 E 428 ~~ 0.054 
818 202 5.07 E 432 0.204 
855 73 2.51 E 450 0.071 
936 Se 3.44 E 494 0.046 
939 87 1.21 E 497 0.077 
976 628 2.69 E 513 0.535 
981 84 5.88 E 517 — 0.071 
999 12 1.14 E 528 0.060 
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Table I 


Root of x? + x? —2x —-1 


21 
= 2cos— 
Oy cos 
BA Ty Ue) 


420 2 3 
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1 11 
3 kee 2 ee 
ZS 12a 
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Sie ta 
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620 
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r.2 4a 


a |= sa | 
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1 
15 


1% 


- 6 8 
tL 


Leto 19" 9 
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| Jig} 3} TUE 
Bl Gl 
*m 4 7 


93 12 


y, 


1 


5 
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i yp 
9 
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2 2 


mies 32h 
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680 
700 
720 
740 
760 
780 
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coe? eee, 1 
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PRT MO RI. 7) 


Dies I 
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1*p 4 11 19 
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24 2 4 3 
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| ah ee 
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Table I (continued) 


_—_—_— eS 


2, 
Obs =2eos =, Root oii. 4 x? — 2x — | 


fo 34 85> 67 8 9 10 P23) 4 S56 88 AG 
800 eo lett t 1 201 f 1 6 «6 6 2d 
820 cee 2 eee LO *g - ft 110 1 3 «1:14 Peas 
840 Sal cee eG oe Ll! 1 OS 9 tr OT I 26 
i) 6 NM Aaa? Si SI, es U2 a a Bate 2 1. Chee Z ees 
880 Sete 3) Sele. 4. Of 3.) 13 1 4 20 ae 
900 S10 >7aer *s 2081 1 1 66 026 L444 6 5 5001 1 
920 oo lle 6 bas 1 lt 1 TD eles. 3: 9) et 
940 aed 1 2, Ay ie 3 Sl ao 3 <4) 2 Sa 
960 io? 14s 41 =66 46 L246 3 TF 3 Rae 7 
980 eee Ss lied 1 64552 be 2 2 1. 72a 
a = 636 b= 119 €=425 ad=202 e = 136 f = 699 
g = 424 h = 165 t= 114 j = 283 k= 267 m = 716 
n= 108 p=704 q=904 r= 124 s= 152 
Table II 
, 2n 
a, = 2cos—, Root of x? + x? — 2x -1 
Frequency Counts 
1 401 15 3 Ceol 59 2 136 «1 
2 168 16 4 341 60 1 52a 
3 109 17 > 36 s6. 1 oa 1 165 1 
4 60 18 3 2 | oe 2 202. —s 1 
5 40 19 Veo 33 | 66 «1 267 «1 
6 23 20 «8 44 1 a 1 283 1 
7 19 pa 3) -45 1 im 3 424 1 
8 is 225e3 46 1 87 «1 425,51 
9 19 vy. a) 49 1 9 «(1 636 «61 
10 12 2503 50 2 a 699 1 
11 10 20ees ene | 108 =i 704 «1 
il94 7 27a eee 114 1 716 et 
13 12 30 «4 mish ed 119 1 904 1 
14 | 52) oe ie ae! 124.1 
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Table III 


A 
Oy =2eos=, Root of x? + x? — 2x —1 


n a,, Qn Tn 
31 94 1.34E 24 1.854 
57 60 2.49 E 29 0.958 
60 50 2.058, 31 0.739 
65 53 5.83E 34 0.698 
78 636 2.51E 42 6.978 
93 119 3.20E 51 1.054 
105 62 196E 58 0.480 
120 425 161E 68 2.815 
243 =e) 7.83 E 124 0.196 
279 202 4.61 E 144 0.617 
286 136 3.97 E 148 0.404 
BPA 62 3.30 E 166 0.166 
357 - 699 1.99 E 185 1.665 
381 o7 1.44 E 200 0.213 
420 61 2.64 E 220 0.121 
463 56 4.85 E 244 0.101 
483 424 . 6.1S E255 0.728 
492 165 6.02 E 260 0.278 
496 114 3.47 E 263 0.191 
514 72 4.03 E 272 ~ O.D16 
583 283 3.59°E 305 0.406 
594 267 1452312 0.375 
630 75 4.20 E 331 0.100 
696 716 2.26 E 366 0.856 
700 108 8.11 E 369 0.129 
737 87 2.72 E 387 0.098 
741 Se) 1 SSE cou 0.066 
762 704 2.68 E 400 0.770 
830 904 5.77 E 434 0.909 
855 124 2.41 E 448 0.121 
903 57 3.80 E 472 0.053 
905 152 5.92 E 474 0.141 
910 66 2.45 E 478 0.061 
918 50 2.34 E 484 0.045 


948 a9 1.04 E 500 0.052 
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Table I 


a, Root of x*—x—1 
6 7 8 9 10 
2) 2h!) Dp) 
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2) Eee 
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420 


440 


te 3 1S 2a 9 


1 
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i i6 
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480 
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2 26 21 


139 4 
3242. 2 
22 eee 3 
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560 
580 


18 6 *f 6 32 oie I ae a ae > 2 


y2 442 


600 
620 


640 


1 


e42 17-1 


2-2 3.3 14 
PAS 2) 


imei 5S lel i 2 2 


3 
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Sezer 


660 
680 
700 


720 


117 4-302 


2) 
4 5 416 3 


ee ers ae 


Pe2ouls. 35 


1 34 
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1 
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Table I (continued) 


a 


a,, Root of x*-x—-1 


12-3: (Soe Sao 7 ae hor 3 4 S56 oP OC 
RR a a ee 
800 13 3 6°1 “74147 2eee 227% Oo 2a es A 
820 £19 2 1 2. ‘Oye eeo 4°11 4 3 25 f 1 "ee 
840 244 3 150 le? ie 05 1 1 *po 7S 2) 1 gee 
860 Sev 1 dod We26ee2 24g Py s.2, Tere? a 2 eee 
880 tS 1 6 | eae 3 1 2 1 4.12 21 (9°52. 1 oe 
900 9°24. 1.393) 53> Seer 1 1 ese 2 Oe ae 
920 Pipe ca) ol Vee O27 3.11435 t 4 1 6 ore 
940 RS Sas or oa See ee te ee | 2) 3. 9 7 -1 732) 8 1 hae 
960 | URES Et came Us ex UR HS Uo ff 301) 79 e207 2 Gee 
980 6 4 ol 16 2s" t 3. 2 8 1310 £ 30 1 15 0tee 


a— 76) b= 195 c = 166 d— 264 e= 7/01 f = 11644 
g = 169 h = 673 b= 457 j = 409 k = 274 m= es 
n= 124 p=172 q = 1033 r=110 s = 684 t = 1292 


Table II 


a, Root of x*—x-1 
Frequency Counts 


1 406 16 6 2) Ra | 48 1 169 «$1 
2 162 17 6 So 3 50. 1 ig 1 
3 89 php j/ S53 3 =e | 174 1 
4 67 1 341 54. 195 1 
5 41 20002 55, | ee 264 = 1 
6 28 2b 4 26 «1 > an 274 = 1 
“i 23 eon lt st 69 «(1 409 I 
8 21 2x. 63 Sho tt 7 1 457 61 
9 11 24 «3 39 «(1 74 #1 673 i! 
10 14 25 4 41 1 85 2 684 1 
11 11 26S 42 2 91 1 701 1 
2 16 4 (| 43 2 99 1 761 1 
13 6 281 44 1 mo 61 1033 1 
14 6 Zoe 46 1 124 1 1202 ~=1 
15 3 30.2 47 1 166 1 11644 1 


[APP. B, 5] COMPUTATIONAL METHOD 121 
Table III 
a Root of x — x — 1 
WW a, Qn ty 
151 ao Zo0 Be 77 0.344 
178 91 3.13 E 91 0.444 
196 5) 1295 102 0.239 
210 761 5.27 E 109 3.096 
POM 55 1.58 E 118 .208 
269 18) 2.00 E 138 0.235 
28 195 3.54 E 141 0.609 
274 54 1.91 E 143 0.166 
285 9 1.35 E 149 0.294 
320 85 8.07 E 163 0.230 
335 85 Bee ee ae 5314) 0.219 
345 166 L387 E177 0.413 
366 264 4.47 E 187 0.621 
382 69 — 1.92 E 195 0.157 
455 701 1.44 E 229 1.347 
588 11644 6.58 E 300 17.042 
643 169 7.49 E 327 0.226 
pe 673 8.41 E 369 0.797 
762 457 1.47 E 391 ORS 
806 409 1.39 E 413 0.434 
807 74 1.03 E 415 0.078 
812 274 4.43 E 421 0.284 
839 174 8.98 E 437 0.174 
840 124 1.11 E 440 0.123 
845 50 2.03 E 444 0.050 
855 172 9.54 E 449 0.167 
870 1033 3.09 E 461 0.979 
897 a) 3.27 E 474 0.048 
908 110 2.20 E 482 0.100 
944 684 3.68 E 503 0.594 
963 1292 UO 1S 1.096 
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Table I 
x® — 9x4 — 4x3 + 27x? — 36x — 23 


a3 = 2+ 3, 
2345 6 7 8 9 
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Table I (continued) 
ay = /2+./3, x®—9x* — 4x3 + 27x? — 36x — 23 


12 38 5 6 7 8 910 2 eS 4a Ao) eso 10 
800 i ets oo oeee lao tf Sati 2 2 12 5 1 7 
820 moe by Seetoeiseetees -§ 1 oper, 3 2 44 1 1 
840 a 1 eerie 2 8 OS 28 1 1 AT a A OB 
860 2047 2. 3° 145 ei 29 1 EOS os 3 al 
880 Socleet. 1. 2 Saige? 2 2 Dot tf 2 43 2 esa 
es eee Geet 2 1 663.0) 6 3 3 06Ud3 23 Pl 1 
920 Sieve 9° 1 2 LS? 1.29 1 Ss 2 2 
940 2 Oot SU eee YL SB 9 tiem es 3. 11. Sale 1 
960 So) as 2 tt 6 1 BU ie 335 eel 2 
980 See eee, 1 8 22 le hl; 25 ey 10 
a= 123 b = 1446 e=126 d=121 e= 154 j =452 
g = 315 h=135 i=4103 j=120 k=331 m = 184 
n= 133 p = 430 q = 298 r= 150 s = 208 t = 186 
u = 138 
Table II 
a3 = 3/24 ./3, x® — 9x4 — 4x? + 27x? — 36x — 23 
Frequency Counts 
1 418 20 1- 42 1 iV) 1 
2 156 21 2 44 1 126 1 
3 105 a 3 47 1 133 i 
4 56 23 4 50 1 135 1 
5 38 24 2 «53 1 138 1 
6 20 25 2 25 1 150 1 
7 30 26 1 58 1 154 1 
es, a, 3 61° 1 184 1 
9 12 28 2 63 1 186 1 
10 11 29 3 69 1 208 1 
11 8 207) 5 74 I 298 1 
12 7 ail 1 fe) 1 315 1 
13 11 az 1 80 1 ast v 
14 4 33 1 85 1 430 1 
15 4 34 2 94 1 452 1 
16 5 a5 1 99 1 1446 1 
17 4 36 2 103 1 
18 4 37 3 120 1 
19 4 38 4 121 1 
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Table ITI 
a, = 3/2 + ./3, x® — 9x* — 4x3 + 27x? — 36x — 23 
n ; ay, Qn tn 
3 i123) (Leet ye A Su 

25 69 145E 14 2.462 

54 1446 3.93 E 28 24.700 

83 74 2IQE 43 0.779 
102 126 1225S 1.074 
121 85 : 136E 63 0.608 
209 121 8.04 E 107 0.497 
229 154 3.52 E 118 0.578 
231 61 6.49 E 120 0.224 
260 452 9.82 E 135 1.473 
295 63 1.06 E 154 0.180 
324 315 6.59 E 169 0.818 
331 135 : 4.31 E 174 0.341 
364 103 4.51 E 189 0.240 
374 94 1.00 E 196 0.212 
581 99 2.74 E 294 0.148 
619 120 6.52 E 313 0.168 
621 331 2.18 E 316 0.459 
698 184 9.56 E 356 0.226 
TAY 75 7.91 E 364 0.091 
739 133 9.57 E 377 0.154 
753 58 3.86 E 386 0.066 
767 430 ; 2.47 E 398 0.472 
772 298 4.13 E 402 0.324 
804 150 1.49 E 418 0.158 
806 5) 2.48 E 420 0.057 
899 50 2.45 E 463 0.048 
911 53 6.77 E 470 0.050 
927 208 5.44 E 480 0.189 
945 80 9.75 E 491 0.071 
969 186 4.09 E 505 0.161 


975 138 2.61 E 509 0.119 


The program to do the calculation was written in Fortran, using 
machine-language subroutines for multiple-precision integer arithmetic to 
handle the coefficients of the polynomials. The calculation of a, was 
done in floating-point arithmetic (approximately 14 significant digits), 
using a floating-point approximation to P,(x) (suitably scaled). This pro- 
cedure avoids the use of multiple-precision arithmetic in any trial-and- 
error steps, and so makes for greater efficiency. One could be even more 
efficient, using an idea suggested by Lehmer [5], and compute several 
successive partial quotients from an approximation to P,(x). It is possible 
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to find P,in(x) from P,(x) and 4@,,4,41,.--:4,+m—1 With less multiple- 
precision calculation than is needed to find all the intervening poly- 
nomials explicitly. The additional complication in the program, however, 
did not seem worthwhile, since the results given here were obtained by 
the simpler method in about 6 minutes on an IBM 360/91. A listing of 
the actual program may be obtained on request from Trotter. 
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APPENDIX C 


Addendum to 
Continued Fractions for 
Some Algebraic Numbers 


By S. LANG at New Haven and H. TROTTER at Princeton 


References [2], [3], and [4] came to our attention after the proof-sheets 
of [1] had been corrected. It is clear that the computational method we 
used is essentially the same as that used in [2] and [4], and is presum- 
ably the same as that used in [3] (which does not give details of the 
computation). 

Reference [4] gives the results of a y?-test comparing the observed 
frequencies of partial quotients of certain algebraic numbers with the 
theoretical frequencies for a “random” number. Results are reported for 
nine algebraic numbers, for each of which between 700 and 800 partial 
quotients were calculated. Nothing was found to suggest non-randomness 
except for a very low value of x? (indicating unusually good agreement 
between expected and observed frequencies) for the expansion of the cube 
root of 2. After some discussion the authors remark “... the impression 
persists that the expansion of 2’? is peculiar. Probably the expansion 
will have to be carried to many more terms to verify or contradict this 
impression.” 

It therefore occurred to us that it might be worthwhile to exhibit 
the results of applying a similar y?-test to the expansions that we had 
calculated. Following [4], we divided the partial quotients into ten 
groups consisting of 1, 2, 3, 4, 5 and 6, 7 and 8, 9 through 12, 13 
through 19, 20 through 40, and over 40. For each of the eight numbers 
for which we obtained expansions, we give the y* value obtained from 
the distribution of the first 1000 partial quotients, and in the column 
headed P, the approximate probability that the y? for a random sample 
would be no larger. (The probabilities are computed for the ordinary 


Reprinted from J. reine angew. Math. 267 (1973), pp. 219-220. 
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z’-distribution on nine degrees of freedom. This is not strictly correct 
because the partial quotients of a “random” number are not independent. 
The error involved is assessed in [4], and we agree with the authors that 
it is negligible for present purposes.) For the first six numbers (the 
numbers of degree 3) we give the same information for the distribution of 
the first 3000 partial quotients. 

The rows of the table correspond to the numbers reported on in [1]. 
Thus the first five are the cube roots of 2, 3, 4, 5, and 7, and the last 
three are the positive roots of 


and 
9 = 4x x? = 36x — 23: 


The results do not suggest any significant departure from random 
behavior. In particular the anomaly observed in [4] for the cube root of 
2 appears not to persist when the expansion is carried further. 
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ics without a very extensive building of theories. One general idea is 
that algebraic numbers will exhibit a behavior that is the same as 
almost all numbers in a probabilistic sense, except under very spe- 
cific structural conditions, namely quadratic numbers. Results for 
almost all numbers (due to Khintchine) show an interplay between 
calculus and number theory, which will also show undergraduates 
how analysis mixes with number theory. 
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